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PREFACE 


This  thesis  Is  concerned  vlch  the  detection  of  carbon  monoxide  In 
the  atmosphere  at  various  altitudes  and  temperatures  from  a  shuttle- 
borne  differential  absorption  Ildar.  It  parallels  an  earlier  study  by 
Robert  H,  Wank,  which  focused  on  carbon  dioxide,  and  Is  Intended  to  be  a 
companion  volume  to  his  study. 
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various  difficulties  encountered  along  the  way.  Also,  a  note  of  thanks 
Is  In  order  for  my  sponsoring  organization,  Foreign  Technology  Division, 
to  whom  I  am  grateful  for  the  opportunity  to  complete  my  degree  program. 
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Abstract 
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\ 

^‘^mote  sensing  of  carbon  monoxide  from  a  shuttle-borne  differential 

absorption  Ildar  (DIAL)  was  Investigated.  The  primary  Ildar  wavelength 

rv\ '  c.  r  rii  ", 

Is  In  the  Infrared  at  4.64yii^  and  takes  advantage  of  a  spectral  coin¬ 
cidence  between  the  R(2)  line  In  the  fundamental  absorption  band  of  CO 
and  the  frequency^oubled  R(18)  emission  line  of  the  002*  laser. 
Extinction  coefficients  for  CO,  H2O,  Rayleigh,  and  Mle  scattering  were 
determined  in  order  to  compute  the  return  signal  strength  from  various 
altitudes.  Direct  detection  was  found  to  be  unsatisfactory,  but  hetero¬ 
dyne  detection  was  found  to  be  suitable  If  shot-averaging  Is  used.  With 
heterodyne  detection  the  system  was  determined  to  be  capable  of  making 


accurate  measurements  of  CO  In  the  troposphere,  but  performance  in  the 
stratosphere  was  found  to  be  marginal^ 


DETECTION  OF  ATMOSPHERIC  CARBON 
MONOXIDE  FROM  A  SHUTTLE-BORNE  LIDAR 


I.  Introduction 


A  shuttle-based  Ildar  system  offers  unique  advantages.  From  Its 
vantage  point  In  orbit ,  the  system  would  have  access  to  a  large  percen¬ 
tage  of  the  earth's  atmosphere  and  be  able  to  make  global  measurements 
with  a  spatial  resolution  and  sensitivity  that  could  be  equalled  only  by 
a  large  network  of  stations  on  the  ground  or  an  extensive  airborne 
program.  In  recognition  of  this  fact,  the  National  Aeronautics  and 
Space  Administration  (NASA)  produced  in  1973  a  report  detailing  shuttle 
payload  requirements  for  AMPS  (Atmospherics,  Magnetospherics,  and 
Plasmas  In  Space),  a  scientific  payload  to  be  flown  in  conjunction  with 
the  European  Spacelab.  Included  in  this  study  was  the  scientific 
rationale  for  a  spaceborne  atmospheric  Ildar. 

In  1977,  the  NASA  Atmospheric  Lldar  Working  Group  was  convened. 

One  of  Che  purposes  of  Che  group  was  to  assist  NASA  in  defining  a 
^  jgram  that  would  eventually  lead  to  the  development  of  the  necessary 
flight  systems  (Ref  1).  Also  in  1977,  the  European  Space  Agency 
released  a  report  detailing  the  technical  aspects  of  a  multidisciplinary 
Ildar  system  suitable  for  shuttle  flight  In  the  Spacelab  payload  (Ref  2) 
This  report  was  used  extensively  by  Kweder  in  an  evaluation  of  remote 
sensing  of  gases  using  Ildar  techniques  (Ref  3),  and  later  by  Wank  in 
his  study  on  the  detection  of  carbon  dioxide  from  a  shuttle-based  Ildar 
(Ref  4).  The  current  work  is  also  based  partially  on  the  results  of  the 


European  Space  Agency  report ,  and  should  be  viewed  as  a  companion  volume 


to  the  studies  by  Kweder  and  Wank. 

This  report  addresses  the  feasibility  of  a  shuttle-borne  Ildar 
system  designed  to  measure  atmospheric  carbon  monoxide.  Carbon  monoxide 
(CO)  Is  a  toxic  pollutant  and  Its  effects  on  urban  areas  where  the  auto¬ 
mobile  density  Is  high  have  been  studied  extensively.  Less  well  known, 
however,  are  Its  effects  on  the  earth's  atmosphere  taken  as  a  whole. 

For  Instance,  CO  Is  known  to  figure  prominently  In  the  production  and 
transport  of  other  pollutants,  such  as  NO,  O3,  and  NO2  (Ref  5:855). 

Also  of  Interest  are  the  atmospheric  sources  and  sinks  of  CO.  Direct 
air  sample  measurements  Indicate  that  the  primary  sources  of  CO  are  the 
combustion  of  fossil  fuels  and  the  ocean,  where  bacterial  decomposition 
of  plankton  Is  suspected  (Ref  6:2217).  The  manner  In  which  CO  Is 
destroyed  Is  not  well  understood.  Oxidation  by  OH,  H2O2,  and  HO2  radi¬ 
cals  In  the  stratosphere  accounts  for  part  of  the  process ,  but  there  Is 
also  a  tropospheric  sink  whose  nature  Is  unknown  (Ref  7:2866). 

There  Is  also  evidence  to  suggest  that  trace  gases  In  the 
atmosphere  may  have  profound  long  term  effects  on  the  earth's  weather 
patterns.  For  example,  an  Increase  In  carbon  dioxide  levels  could  lead 
to  a  dramatic  Increase  In  the  average  global  temperature  because  of  what 
has  been  termed  the  "greenhouse"  effect .  The  manner  In  which  CO  might 
effect  changes  In  weather  patterns  Is  not  clear  and  more  data  are 
necessary  If  accurate  predictions  are  to  be  made  of  the  long  term 
effects  of  the  burning  of  fossil  fuels  as  a  primary  source  of  energy. 
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earth's  atmosphere,  global  measurements  of  the  CO  number  density  as  a 
function  of  altitude  above  the  earth's  surface  are  necessary.  A 
shuttle-based  Ildar  system  Is  Ideally  suited  for  this  task. 


Also  addressed  In  this  study  Is  the  feasibility  of  a  shuttle-borne 
Ildar  system  to  detect  and  measure  hot  CO  concentrations  near  the  sur¬ 
face  of  the  earth.  Such  a  capability  could  be  used  to  analyze  the 
emissions  from  factory  stacks  or  perhaps  to  locate  and  characterize  the 
exhaust  plumes  from  gas  turbine  or  rocket  engines. 

Problem 

The  problem  of  this  thesis  Is  to  Investigate  the  feasibility  of  a 
shuttle-borne  Ildar  system  designed  to  detect  and  accurately  measure 
atmospheric  carbon  monoxide  concentrations  at  all  altitudes.  The  basic 
parameters  of  the  system  are  those  recommended  in  the  European  Space 
Agency  study  (Ref  2).  The  primary  wavelength  selected  for  the  system 
utilizes  a  laser/absorption  line  coincidence  reported  by  Killlnger  et  al 
(Ref  8:403). 

General  Approach 

The  Investigation  Is  Initially  driven  by  the  choice  of  primary 
wavelength.  Possible  candidate  Ildar  systems  that  can  be  used  to  pro¬ 
vide  the  required  data  are  analyzed,  and  a  primary  candidate  Is 
selected.  The  number  of  variables  In  the  analysis  Is  reduced  by 
assuming  that  the  basic  telescope  and  associated  optics  are  fixed.  The 
system  variables  then  are  those  of  the  type  of  laser  used  and  the 


technique  used  to  detect  and  process  the  received  signal.  A  final 


system  configuration  Is  selected  on  the  basis  of  the  analysis  results 
and  Its  performance  Is  assessed. 

Assumptions 

The  atmospheric  carbon  monoxide  Is  assumed  to  be  uniformly  mixed  at 
0.075  ppmv  (parts  per  million  by  volume),  that  is  this  relative  con¬ 
centration  at  sea-level  Is  the  same  at  all  altitudes  up  to  30  km. 

The  atmosphere  Is  considered  to  be  nonturbulent  and  without  cloud 
cover.  The  temperatures,  pressures,  and  total  molecular  number  den¬ 
sities  are  those  given  In  the  U.S.  Standard  Atmosphere,  1976  (Ref  9). 
The  water-vapor  profile  Is  also  taken  from  this  reference. 

Aerosol  concentrations  are  from  the  model  by  McClatchey  (Ref  10). 
Visibility  Is  assumed  to  be  25  km  at  0.55  microns. 


II.  Lldar  Detection  of  Carbon  Honoxlde 


Lldar  techniques  have  been  used  to  monitor  the  earth's  atmosphere 
since  the  early  1960's.  The  technology  has  progressed  to  the  point 
where  high-performance  systems  can  now  be  contained  In  a  relatively  com¬ 
pact,  light-weight  structure  suitable  for  use  on  a  space  vcehlcle.  The 
purpose  of  this  study  Is  to  assess  the  feasibility  of  one  such  system, 
designed  to  measure  atmospheric  carbon  monoxide  from  the  payload  bay  of 
the  space  shuttle. 

The  acronym  "lldar"  stands  for  light  Retention  and  ranging.  Lldar 
systems  are  also  sometimes  called  laser  radars.  They  function  very  much 
like  their  microwave  counterparts :  electromagnetic  radiation  Is 
transmitted  Into  the  atmosphere  by  a  pulsed  or  continuous -wave  (CW) 
laser  source,  the  laser  energy  Is  absorbed  and  scattered  by  gas  molecu¬ 
les  and  aerosols  In  Its  path,  and  a  portion  of  the  scattered  radiation 
Is  then  detected  and  measured.  The  amplitude  and  other  characteristics 
of  the  scattered  radiation  that  enters  the  receiver  can  be  used  to  Infer 
various  facts  about  the  volume  of  atmosphere  under  consideration.  For 
Instance,  lldar  systems  can  determine  absorption,  aerosol  content,  and 
also  number  density  of  the  various  constituents  along  the  laser  path. 

Various  lldar  techniques  and  system  configurations  are  used  to  make 
atmospheric  measurements.  If  the  laser  transmitter  Is  colocated  with 
the  receiver,  the  system  Is  referred  to  as  single-ended.  If  the 
transmitter  and  receiver  are  separated,  the  system  Is  termed  double- 
ended.  The  techniques  currently  In  use  can  be  separated  Into  three  main 
groups:  Raman,  resonance  emission  (fluorescence),  and  differential 


absorption  (DIAL).  As  one  might  expect,  Raman  Ildar  attempts  to  detect 
molecular  Raman  frequency  shifts  In  the  returned  signal.  Each  gas  mole¬ 
cule  has  Its  ovm  unique  Raman  spectrum.  Therefore,  In  principle,  such  a 
system  could  use  a  single  exltatlon  frequency  to  obtain  Information 
about  all  the  different  gas  molecules  along  the  laser  path.  The  biggest 
disadvantage  to  Raman  Ildar  Is  the  extremely  small  scattering  cross- 
section  for  the  molecules  of  Interest  (typically  10“^^  cm2  sr"!).  This 
limits  the  useful  range  of  the  Raman  technique  to  a  few  km  (Ref  11: 
1644). 

Resonance  emission,  or  fluorescence,  Ildar  techniques  make  use  of  a 
tunable  laser  source  to  excite  selectively  the  various  molecules  of 
Interest.  An  Infrared  laser  can  be  used  to  excite  vibrational- 
rotational  transitions  of  the  molecule,  or  a  laser  operating  In  the  UV 
or  visible  can  be  used  to  directly  excite  atomic  transitions.  The  mole¬ 
cule  or  atom  subsequently  fluoresces  and  the  resulting  radiation  Is 
detected  and  analyzed.  The  technique  Is  simple  and  more  sensitive  than 
Raman,  since  the  required  cross-sections  are  much  larger.  The  only  real 
problem  with  this  technique  In  the  Infrared  Is  what  Is  known  as 
quenching.  That  Is,  non-radlatlve  decay  processes  such  as  frequent 
collisions  occur  more  quickly  than  radiative  decays.  As  a  result,  the 
number  of  radiative  emissions  becomes  uncertain  and  the  accuracy  of  the 
measurement  Is  called  Into  question  (Ref  4:15). 

The  last  technique  to  be  discussed  Is  differential  absorption 
Ildar,  or  DIAL.  In  this  technique,  radiation  Is  transmitted  Into  the 
atmosphere  at  two  distinct  wavelengths:  one  of  which  Is  tuned  to  a  par¬ 
ticular  absorption  line  of  the  molecule  of  Interest,  while  the  other  Is 
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slightly  offset  from  the  absorption  line  to  provide  a  reference  signal. 
The  backscattered  return  signals  at  the  two  wavelengths  are  then  com¬ 
pared  and  analyzed.  The  amplitude  difference  between  the  two  Is  a 
direct  Indication  of  how  much  laser  energy  was  absorbed  by  the  gas  mole¬ 
cules  of  Interest,  and  hence,  their  number  density  can  be  determined. 

If  a  pulsed  laser  transmitter  Is  used,  then  the  resulting  return  signals 
can  be  compared  as  a  function  of  time,  and  range-resolved  measurement 
data  then  become  available.  The  differential  absorption  technique  Is 
very  sensitive  and  can  be  used  over  long  ranges,  provided  sufficient 
laser  power  Is  available. 

Of  the  three  Ildar  techniques  Just  discussed,  only  resonance 
emission  and  differential  absorption  appear  suitable  for  use  In  this 
study.  The  resonance  technique  suffers  from  the  quenching  problem  at  IR 
wavelengths,  and  therefore  should  be  used  only  at  very  high  altitudes 
where  the  collision  probability  Is  lower.  Thus,  differential  absorption 
Ildar  appears  to  be  the  only  technique  available  to  measure  the 
atmospheric  content  of  trace  gases  such,  as  carbon  monoxide,  over  the 
broad  range  of  altitudes  required  for  a  shuttle-based  system  (Ref  12: 
624). 


System  Definition 


The  Ildar  system  that  will  be  considered  in  this  study  Is  single- 
ended  and  utilizes  the  differential  absorption  measurement  technique. 

The  primary  wavelength  that  will  be  used  Is  recommended  by  Kllllnger  et 
al  (Ref  8:403).  The  wavelength  is  at  4.64^  in  the  Infrared  and  utilizes 
a  spectral  coincidence  between  the  R(2)  transition  In  the  fundamental 


absorption  band  of  the  CO  molecule  at  2154.5960  cm~*,  and  the  frequency- 
doubled  R(18)  laser  line  in  the  00*1-02*0  band  of  CO2  at  2154.6050  cm”^. 
The  reference  wavelength  must  be  outside  the  R(2)  line  of  CO,  but  close 
enough  to  the  primary  so  that  the  difference  in  absorption  experienced 
at  the  two  wavelengths  can  be  directly  attributed  to  the  R(2}  line.  The 
wavelength  selected  for  this  study  uses  the  frequency -doubled  R(10)  line 
of  the  CO2  laser  at  2143.7674  cm"^  (Ref  13). 

The  spectral  coincidence  involving  a  frequency -doubled  transition 
of  the  CO2  laser  is  fortunate.  Carbon  dioxide  laser  transitions  are  very 
accurately  known  and  CO2  laser  system  technology  is  well  advanced.  In 
addition,  their  efficiency  is  typically  10-20  percent,  which  is  well 
above  average. 

For  a  pulsed,  high  peak  power  application  such  as  that  required  in 
this  study,  a  transversely  excited  atmospheric  pressure  (TEA)  CO2  laser 
should  be  used.  Because  of  the  transverse  excitation  scheme,  TEA  lasers 
automatically  emit  giant,  Q-switched  like,  pulses  due  to  a  gain 
switching  action  in  the  lasing  medium.  The  pulses  are  nominally  100  ns 
long,  and  their  repetition  frequency  can  be  controlled  by  the  excitation 
process.  Pulse  energies  of  up  to  20  joules  per  100  ns  pulse,  at  a  repe¬ 
tition  rate  of  up  to  1200  Hz  are  reported  (Ref  14:110). 

Carbon  dioxide  TEA  lasers  have  other  advantages  as  well.  They 
operate  at  room  temperature  and  usually  at  atmospheric  pressure.  They 
are  also  self-contained  and  do  not  require  venting  of  exhaust  gases. 
Finally,  they  are  line -tunable  over  a  range  of  about  80  or  so 
rotational-vibrational  transitions . 


They  do  have  some  problem  areas  that  should  be  noted.  The  short 
duration,  high  energy  pulses  will  cause  electromagnetic  interference 
(EMI)  problems  on  the  shuttle,  if  not  properly  shielded.  The  normal 
operating  voltage  for  a  typical  Installation  is  about  40  kV.  This  can 
be  extremely  dangerous  to  operating  personnel.  The  last  problem  area  is 
one  that  will  affect  the  performance  of  the  system.  The  output  spectral 
llnewldth  of  the  CO2  TEA  laser  is  about  2-3  GHz  (Ref  2:78  and  Ref 
14:110).  This  is  about  the  same  order  of  magnitude  as  the  width  of  a 
typical  CO  absorption  line  at  sea-level.  It  is  disadvantageous  because 
it  implies  that  a  detailed  absorption  cross-section  as  a  function  of 
wavelength  must  be  known  if  accurate  measurements  of  number  density  are 
to  be  extracted  from  the  data.  This  point  will  be  fully  addressed  later 
in  the  study.  For  now,  it  is  convenient  to  state  that  the  CO2  TEA  laser 
output  llnewldth  assumed  in  this  study  is  3  GHz  (0.10  cm"^). 

In  order  to  achieve  the  output  wavelength  specified  for  the  system, 
the  emissions  from  the  primary  and  reference  CO2  TEA  lasers  must  be  fre¬ 
quency-doubled.  This  can  be  accomplished  through  the  use  of  a  crystal 
of  CdGeAs2.  Klldal  has  irradiated  this  material  with  150  ns  CO2  laser 
pulses  and  obtained  200  mJ  of  doubled  radiation  with  a  conversion  effi¬ 
ciency  of  lOZ.  Hlnkley  states  that  this  crystal  should  be  capable  of 
providing  several  Watts  of  average  power  and  hundreds  of  mlllijoules  of 
pulsed  power  (Ref  15:251). 

The  transmit /receive  optics  that  will  be  used  in  this  study  are 
those  recommended  in  the  European  Space  Agency  study  (Ref  2:28).  A 


V.'-* 

--.V/.J 


.  m  V 

•f/.'CN*!' 


single  telescope  will  be  used  for  both  functions.  The  telescope  has  a 
diameter  of  100  cm,  and  therefore  a  transmit /receive  aperture  area  of 
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0.7854  m^.  It  has  a  radian  field-of-vlew  (FOV)  of  10"^,  which  equates 
to  a  solid-angle  FOV  of  7.854  x  10“^  sr.  The  efficiency  of  the  receive 
optics  is  assumed  to  be  0.3,  which  includes  all  optical  components  prior 
to  the  detector. 

A  block  diagram  of  the  system  is  shown  in  Figure  1,  and  the  basic 
parameters  selected  for  the  study  are  summarized  in  Table  I  below. 


Table  I 

Lldar  System  Parameters 


Primary  wavelength 

2154.5060  cm-1 

Reference  wavelength 

2143.7674  cm~l 

Laser  pulsewidth 

100  ns 

laser  spectral  llnewldth 

3  GHz 

Laser  overall  efficiency 

ZOZ 

Diameter  of  optics 

100  cm 

Area  of  optics 

0.7854  m2 

Efficiency  of  receive  optics 

0.3 

Efficiency  of  doubling  crystal 

0.1 

Radian  FOV  (full  angle) 

10"^  radians 

Steradlan  FOV 

7.854  X  10~9  sr 

ecrccroi* 


Block  Diagram 


The  Lldar  Equation 


The  primary  equation  that  will  be  used  to  evaluate  system  perfor¬ 
mance  Is  known  as  the  “Ildar  equatlon“.  It  relates  the  received  signal 
power  from  a  scattering  cell  to  that  Initially  transniitted  by  the  laser. 
For  the  single-ended  DIAL  system  used  in  this  study.  It  may  be  written 
in  the  following  form  (Ref  15:76): 


Pi  (/t) 


€  A 


(/) 


where 


Is  the  received  power  from  the  scatter  cell  at 
range  R,  In  Watts 

4  has  value  I  or  2,  and  refers  to  the  primary  and 

reference  wavelengths,  respectively. 

Is  the  transmitted  laser  power  In  Watts 

6  Is  the  efficiency  of  the  doubling  crystal 

K  Is  the  efficiency  of  the  receiver  optics 

7«  Is  the  receiver  integration  time  In  seconds 

c  Is  the  speed  of  light  In  km/sec 

Is  the  atmospheric  volumetric  backscatterlng 
coefficient  In  km~Ur"^ 


A. 


Is  the  effective  area  of  the  receive  optics  In  m‘ 


R  Is  the  distance  from  the  Ildar  to  the  scatter 

cell  In  m 


Is  the  two-way  atmospheric  transmittance  along 
the  path  from  the  Ildar  to  the  scatter  cell 


The  coordinate  system  that  vrlll  be  used  In  this  study  has  Its  ori¬ 
gin  at  the  Ildar  system  on  board  the  space  shuttle,  which  Is  assumed  to 
be  In  a  200  km  circular  orbit.  The  system  Is  pointed  directly  downward 
so  that  the  range  Is  measured  along  a  line  from  the  shuttle  through  the 
earth's  center.  Thus,  sea-level  on  the  earth's  surface  corresponds  to  a 
Ildar  system  range  of  200  km,  while  the  range  to  a  scatter  cell  at  an 
altitude  of  30  km  above  the  surface  would  be  170  km.  In  this  study  the 
results  will  generally  be  presented  as  a  function  of  altitude  above  the 
earth's  surface,  rather  than  as  a  function  of  range  from  the  Ildar. 

The  term  (c7^/z)  In  equation  (1)  has  dimensions  of  length  and 
defines  the  range  resolution  of  the  system,  /i/i  .  During  the  receiver 
Integration  time  7"^  ,  the  laser  pulse  Illuminates  a  geometrical  length 
c  of  the  atmosphere.  However,  since  the  received  signal  must  tra¬ 

vel  a  two-way  path,  the  atmospheric  length  from  which  signals  are 
received  during  this  time  Is  just  half  this  value.  In  this  study 
Is  taken  to  be  1  km,  which  corresponds  to  a  receiver  Integration  time  of 
6.67  microseconds. 

The  volumetric  backscatter  coefficient  In  equation  (1)  has  dimen¬ 
sions  km~^  sr~^.  This  term  quantifies  the  fraction  of  transmitted  power 

that  Is  scattered  back  to  the  receiver  per  unit  scatter-cell  length  per 

-* 

unit  solid  angle.  Here  the  solid  angle  Is  given  by  the  fs  term 
which  Is  now  recognized  as  the  solid  angle  subtended  by  the  detection 
system  from  the  scatter  cell  at  range  R. 

The  Ildar  equation  will  be  more  convenient  to  work  with  If  It  Is 
written  In  a  form  that  Isolates  any  term  that  Is  highly  sensitive  to 
wavelength.  To  this  end,  the  atmospheric  transmission  term  In  equation 


13 


(1)  can  be  written  as: 


(2  ) 


Is  the  two-way  atmospheric  transmission  due  to 
the  R(2)  absorption  line  of  CO 

Is  the  two-way  atmospheric  transmission  due  to 
all  other  factors 


where 


Then,  the  received  power  at  the  primary  wavelength  can  be  written: 


P,  (ft)  = 


P.  € 


r.;/ 


(3) 


and  the  received  power  at  the  reference  wavelength  becomes : 

K  (V) 

where  it  has  been  assumed  that  the  atmospheric  transmittance  due  to  all 
other  factors  except  the  R(2)  line  of  CO  Is  the  same  at  both  wave¬ 
lengths.  This  is  reasonable,  since  the  reference  wavelength  was  chosen 
specifically  to  make  this  assumption  valid. 

It  can  also  be  assumed  that  the  atmospheric  backscatter  coef¬ 
ficients  at  both  wavelengths  are  approximately  equal  to  one  another, 
since  >fl  is  a  slowly  varying  function  of  wavelength  (Ref  2:58). 

Then,  if  the  laser  output  power  is  the  same  for  both  the  primary 
and  reference  systems.  It  can  be  seen  that  the  returned  power  levels 


for  the  two  will  be  related  by: 


P,(li)  =  T^ll  CD 

A 

This  equation  just  expresses  the  fact  that  Px  suffers  selective  absorp¬ 
tion  due  to  the  R(2)  line  of  CO  while  P2  does  not. 

In  a  range-resolved  DIAL  system,  the  above  relation  can  be  used  to 
obtain  an  equation  for  the  average  number  of  CO  molecules  In  a  range 
cell  ,  extending  from  a  range  Rj  to  R2.  The  ratios  of  power 
received  at  the  primary  wavelength  to  that  at  the  reference  at  Rj  and  R2 
are : 

p.m.) 

Pt(p,) 

and: 

P.C/f^) 

Now,  the  atmospheric  transmittance  can  be  expressed,  using 
Bouguer's  (Beer's)  Law,  as  (Ref  15:49): 

a 

=  ««/>  {-i  f 

where  Is  the  CO  number  density 

7"^^  Is  the  absorption  cross-section  due  to  the  R(2) 
line  of  CO  at  the  primary  wavelength 


- 

- 


Then,  Ic  can  be  shown  that: 
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or : 


C/tf) 


Now  If  the  number  density  and  cross-section  of  CO  are  assumed  constant 
and  equal  to  their  average  values  over  the  range  cell,  then: 


P,  C^x)  Fg  (P, ) 

P^  (r7)  F,  CF,) 


«x/> 
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and  the  average  number  density  can  be  found  from: 
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where  AR  =*  R2“Rl 

Since  the  pulsed  Ildar  return  signal  Is  range  resolved,  the  above 
relation  can  be  used  to  calculate  the  CO  number  density  for  each  range 


cell  of  interest 


III.  Molecular  Structure  of  Carbon  Monoxide 


The  evaluation  of  the  Ildar  equation  requires  a  precise  knowledge 
of  atmospheric  transmission  at  the  wavelength  of  Interest.  The  purpose 
of  this  and  the  following  section  Is  to  lay  the  groundwork  for  these 
calculations,  as  well  as  gain  a  firm  understanding  of  the  manner  In 
which  the  CO  molecule  Interacts  with  the  transmitted  laser  radiation. 


The  Carbon  Monoxide  Molecule 


Carbon  monoxide  Is  a  heteronuclear  diatomic  molecule  whose  charac¬ 
teristics  are  well  known.  Atmospheric  CO  contains  four  principle  Isoto- 


,  ,  12^16-  13-16„  12-18-  .  12-17-  ™  .  v  j  ►  4- 

pic  species :  CO,  CO,  CO,  and  C  0.  The  most  abundant  Is 


12^16, 


CO,  which  comprises  98.652  percent  of  atmospheric  CO  (Ref  16:6). 
Table  11  lists  the  molecular  constants  for  this  particular  Isotope. 


1 2  16 

since  It  Is  principally  CO  that  will  be  detected,  the  remainder  of 


the  study  assumes  that  the  term  "CO"  refers  to  this  Isotope • 


Table  II  ^2  16 

Molecular  Constants  for  C  0 


2169.8196 

13.2918 

y* 

1.1229  X  10-2 

tJe  Zt 

-4.7  X  10-7 

Be 

1.931265 

Oe 

6.1025  X  10-6 

All  values  in  cm“^  (Ref  17:327) 
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This  study  Is  concerned  with  the  CO  molecule's  Interaction  with 
electromagnetic  radiation  In  the  Infrared  region  of  the  spectrum.  The 


Interaction  Is  one  between  the  Incident  photon  and  the  molecule's 
changing  dipole  moment.  Heteronuclear  molecules  have  a  permanent  dipole 
moment  associated  with  them.  Rotation  or  vibration  of  the  molecule 
causes  shifts  In  the  dipole  moment  which  allow  the  molecule  to  absorb 
radiation.  Classical  electrodynamics  predicts  that  such  an  Interaction 
can  occur  for  any  frequency  that  matches  that  of  the  rotating  or 
vibrating  molecular  system.  However,  the  quantum  theory  gives  the 
correct  result  that  such  an  Interaction  can  occur  only  for  certain 
discrete  frequencies  associated  with  transitions  between  the  various 
energy  levels  of  the  system.  These  energy  levels  are  characterized  by 
the  rotational  quantum  number  J,  and  the  vibrational  quantum  number  v. 

Distribution  of  CO  Molecules  In  Vibrational  States 

The  energy  difference  between  a  molecule  In  the  first  vibrational 
state  and  one  In  the  ground  vibrational  state  Is  given  approximately  by 
(Ref  18:95): 

(/J) 

where  the  needed  values  can  be  obtained  from  Table  II,  resulting  In  an 
energy  difference  of  about  2143.2  cm~^. 

This  energy  difference  can  be  used  to  calculate  the  ratio  of  mole¬ 
cules  In  the  first  vibrational  state  to  those  In  the  ground  vibrational 


state.  According  to  the  Maxwell -Boltzmann  distribution  law,  the  ratio 
is  given  by: 


where  k  is  the  Boltzmann  constant  and  T  is  the  absolute  temperature. 

For  a  sea-level  temperature  of  296*K,  the  ratio  is  about  3.0  x  10”^. 
Therefore,  to  a  good  approximation,  all  of  the  CO  molecules  are  in  the 
ground  vibrational  state  at  this  temperature.  In  addition,  the  mean 
atmospheric  temperatures  at  the  altitudes  considered  in  this  report  are 
equal  to  or  less  than  296*K,  so  the  above  statement  can  be  extended  to 
encompass  all  atmospheric  CO  molecules  as  well. 

"Hot"  CO  in  the  atmosphere,  such  as  that  emanating  from  a  factory 
stack  or  a  rocket  engine  exhaust  plume,  might  have  molecular  tem¬ 
peratures  ranging  from  400*K  to  1000*K.  As  an  upper  limit,  consider  the 
N,  / ratio  at  lOOO’K: 

-  z 

¥.S8  *  lo  (/C) 

Even  at  high  temperatures,  most  of  the  molecules  are  in  the  ground 
vibrational  state. 

Absorption  Line  Positions  for  the  1-0  Band  of  CO 

The  fundamental  absorption  band  for  CO  is  characterized  by  tran¬ 
sitions  where  the  initial  vibrational  state  is  the  ground  state  v»0, 
and  the  final  state  is  the  v»l  state.  The  transitions  are  further 


characterized  by  the  value  of  the  rotational  quantum  number  J. 

To  a  first  approximation,  the  molecular  energy  levels  associated 
with  the  transitions  can  be  represented  by  two  additive  components  asso¬ 
ciated  with  (1)  the  rotation  of  the  molecule,  and  (2)  Its  vibrational 
state.  A  very  careful  study  of  the  vibrational -rotational  Interaction 
was  made  by  J.L.  Dunham.  He  expressed  the  energy  levels  of  the  various 
states  as  a  double  power  series  In  J  and  v  (Ref  18:109): 

T C>,J)  =  ^  c,t) 

where  the  y^j  coefficients  in  the  series  are  known  as  the  Dunham  coef¬ 
ficients.  This  series  Is  similar  to  the  usual  spectroscopic  represen¬ 
tation  in  that  5*  ,  etc. 

This  study  will  use  the  Dunham  coefficients  for  CO  tabulated  by 
Todd  et  al  (Ref  19:224)  to  calculate  the  absorption  line  positions  for 
the  1-0  band  of  the  CO  molecule.  The  coefficients  used  In  the  energy 
level  calculations  are  given  In  Table  III.  The  line  positions  calculated 
from  these  energy  levels  are  listed  In  Appendix  A.  These  positions  are 
stated  in  the  Todd  study  to  be  accurate  to  -  0.0005  cm“^.  The  P,R  branch 
terminology  is  standard  and  refers  to  a  »  -1  transition  for  the  P 
branch,  and  a  AT  »  1  transition  for  the  R  branch. 


Calculation  of  Line  Strength  for  the  CO  1-0 


The  line  strength,  or  Integrated  intensity,  for  an  Individual  line 


in  the  CO  absorption  spectrum  Is  defined  according  to  the  following 


2169.813844 


-13.2884214 
10.52536  X  10"3 
5.6644  X  10“^ 
10.0327  X  10"7 
-3.18476  X  10"® 
1.931280880 
1.75043079  x  10"2 
5.6767  X  10"7 
2.3691  X  10"® 
-6.120712  X  10"® 
10.6399  X  10"10 
-1.78578  X  10"10 
5.6676  X  10"12 
-1.62587  X  10-13 


relation  (Ref  15:240): 
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where  Is  the  line  strength  In  ca^  mol“l  cm“l,  and  TCr)  Is  the 

absorption  cross-section  of  the  molecule  at  wavenumber  V  . 

The  Intensity  variation  of  line  strength  as  a  function  of  the  rota¬ 
tional  quantum  number  J  depends  on  the  relative  fraction  of  the  total 
number  of  molecules  In  that  particular  J  state.  This  suggests  that  the 
Individual  line  strength  can  be  calculated  from: 


s 


(/») 


where  la  the  fractional  number  of  molecules  in  state  J,  and  Z 

la  the  sum  of  all  line  strengths  in  the  band,  or  the  "band  sum".  The 
band  sum  for  the  CO  l-O  band  is  reported  to  be  981.3  x  10“20  cm^  mol  ^ 
cm“^,  at  an  absolute  temperature  of  296*K  (Ref  20:793).  The  ratio 

Is  given  by  (Ref  18:125): 
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where  Qj-  ±s  the  rotational  partition  function ;  and  B  is  the  ground  state 
rotational  constant,  which  has  the  value  1.9225  cm”^  (Ref  19:207).  The 
rotational  partition  function  takes  the  form: 
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For  small  values  of  B,  may  be  approximated  by  the  integral  (Bef 
18:125); 
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The  expression  for  then  becomes: 
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In  equation  (23)  It  Is  assumed  that  the  transition  probabilities 
for  the  band  are  the  same  for  all  lines.  Actually,  the  probability  Is 
weakly  dependent  on  J,  In  that  a  particular  line  Intensity  depends  on 
the  average  value  of  J  for  the  upper  and  lower  states  Involved  In  the 
transition.  Taking  this  Into  account,  the  equation  becomes  (Bef  18:126) 
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where  J'  refers  to  the  upper  state  and  J"  to  the  lower. 

The  equation  for  the  line  strength  can  now  be  written: 


(zr) 


where  QvtOr  vibrational  and  rotational  partition  functions, 

respectively,  S,,ct,)1s  the  line  strength  at  296*K,  and  E“  Is  the  energy 
of  the  lower  state  of  the  transition  In  cm'^.  For  CO,  the  vibrational 
partition  function  Is  Independent  of  temperature  (for  the  range  con¬ 
sidered  In  this  study),  and  the  ratio  Q^(Ti>/av(T)ls  equal  to  one  (Bef 
16:4).  Using  the  approximation  for  given  In  equation  (22),  the  ratio 
Q^Cii)/(J,.<T)  becomes  TJ/T,  Therefore,  equation  (26)  can  now  be  writ¬ 
ten  as : 
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IV.  Atmospheric  Transmission 


Solution  of  the  Ildar  equation  requires  accurate  knowledge  of 
atmospheric  transmission  at  the  wavelength  of  Interest,  In  this  case  the 
Ildar  primary  wavelength  at  2154.6050  cm~^.  Recall  that  the  atmospheric 
transmission  term  In  the  Ildar  equation  was  separated  Into  two  factors. 
One  was  the  transmission  due  solely  to  the  R(2)  absorption  line,  which 
can  be  written  (from  equation  (8)): 


It  Is  clear  from  this  equation  that  In  order  to  calculate  the 
transmission,  the  number  density  of  CO  and  the  absorption  cross-section 
of  the  R(2)  line  at  the  primary  wavelength  must  be  known  as  a  function 
of  range.  The  approach  taken  In  this  study  Is  to  break  the  atmospheric 
path  into  R  range-cells,  each  of  length  AR  »  1  km,  which  Is  the  ver¬ 
tical  resolution  of  the  system.  The  number  density  and  absorption 
cross-section  are  then  assumed  constant  over  the  range-cell.  The 
transmission  can  then  be  calculated  for  each  range-cell  and  the  total 
transmission  obtained  from  a  product  of  the  constituent  range-cell 
transmissions  In  the  path  under  consideration.  Therefore: 

R 

T-J/  =  7T 


where  7;o/,  Is  Che  total  two-way  atmospheric  transmission 

due  Co  Che  R(2)  line  of  CO  to  a  range  R 

7;*  )  is  Che  two-way  atmospheric  transmission  along 

a  1  km  vertical  path  at  a  range  Ri 


7’ea  can  be  written  more  explicitly  as: 


where  the  number  density  and  cross-section  are  understood  to  be  eva¬ 
luated  in  the  range-cell  of  interest.  The  atmospheric  transmission  data 
can  Chen  be  presented  in  tabular  form  for  each  range-cell  or,  by  using 
the  conversion  noted  earlier,  for  each  atmospheric  layer.  In  this 
study,  it  is  assumed  that  any  atmospheric  absorption  effects  above  an 
altitude  of  30  km  can  be  neglected.  Therefore,  the  transmission  must  be 
calculated  for  each  of  30  range-cells  extending  from  sea-level  to  an 
altitude  of  30  km. 

The  second  atmospheric  transmission  term  in  the  Ildar  equation 


species  except  the  R(2)  line  of  CO,  Using  Beer’s  law,  this  term  can  be 
written : 
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where  £C»*)  is  the  extinction  coefficient  for  all  atmospheric  species 
except  the  R(2)  line  of  CO.  More  explicitly,  £<<*)  can  be  written: 


where  is  the  extinction  coefficient  due  to  Rayleigh 

scattering 

Is  the  extinction  coefficient  due  to  Mle 
(aerosol)  scattering 

Is  the  molecular  extinction  coefficient  due 

C 

to  CO2  absorption 

Q  Is  the  molecular  extinction  coefficient  due 

*  to  N2O  absorption 

Is  the  molecular  extinction  coefficient  due 

*  to  ozone  absorption 

Q  Is  the  molecular  extinction  coefficient  due 

to  water-vapor  absorption 

Is  the  molecular  extinction  coefficient  due 
to  nearby  CO  absorption  lines 

The  five  molecules  listed  are  the  primary  constituents  of  the 
atmosphere  thought  to  be  responsible  for  attenuation  of  electromagnetic 
radiation  in  the  wavelength  band  considered  In  this  report  (Ref  15:36). 
Not  all  of  them  have  a  significant  Impact  on  the  calculation  at 
2154.6050  cm“^,  as  will  be  seen. 

Accordingly,  ^ can  also  be  written  as  a  product  of  transmissions 
over  Individual  range-cells.  Therefore,  let; 

7-;/  =  tT  ^ 


where 


7;  (Ri )  is  the  two-way  atmospheric  transmission  along 
a  1  km  vertical  path  at  range 


2 

dR)  can  be  written  more  explicitly  as: 

r/<*)  .  <“ 

where  the  extinction  coefficients  are  understood  to  be  evaluated  in  the 
range -cell  under  consideration. 

The  remainder  of  this  section  will  develop  the  theory  and  calcula¬ 
tions  necessary  to  evaluate  equations  (29)  and  (33),  which  can  then  be 
used  to  calculate  the  atmospheric  transmission  for  the  Ildar  equation 
over  any  path  of  Interest. 


Atmosphere  Model 

The  earth’s  atmosphere  is  generally  regarded  as  being  composed  of 
layers,  each  of  which  is  characterized  by  its  temperature  profile.  The 
lowest  layer  is  characterized  by  a  temperature  that  decreases  linearly 
from  sea-level  to  an  altitude  of  about  12  km.  Above  this  altitude  the 
temperature  starts  to  increase  once  again.  This  point  is  called  the 
tropopause.  It  divides  the  atmosphere  into  two  regions:  the  atmosphere 
below  the  tropopause,  called  the  troposphere;  and  that  above  the  tropo¬ 
pause,  called  the  stratosphere. 

The  range  of  altitudes  that  will  be  considered  in  this  study  is 
from  sea-level  to  30  km.  Above  this  altitude  the  particle  and  molecular 
number  densities  are  so  low  that  any  absorption  that  occurs  is 


negligible 


The  basic  atmosphere  model  used  in  this  study  is  from  the  1976  ver¬ 
sion  of  the  U.S.  Standard  Atmosphere  (Ref  9).  The  temperature, 
pressure,  and  total  number  density  profiles  extracted  from  this 
reference  are  reproduced  as  Appendix  B. 

Carbon  monoxide  is  assumed  to  be  a  uniformly  mixed  gas  throughout 
the  altitude  range  considered  in  this  study  (Ref  16:2).  The  mixing  ratio 
is  variously  reported,  depending  on  time  and  place,  as  between  0,05  and 
0.25  ppmv  (Ref  21:1-2).  In  this  study,  a  mixing  ratio  of  0.075  ppmv  is 
assumed,  which  should  produce  conservative  results  (Ref  16:2),  Using 
this  mixing  ratio,  the  total  atmospheric  molecular  number  density,  and 
the  fact  that  0.98652  of  atmospheric  CO  is  the  Isotope  ^2cl6o;  the 
number  densities  for  this  Isotope  can  be  calculated  as  a  function  of 
altitude.  The  results  are  tabulated  in  Appendix  B. 

The  water-vapor  profile  assumed  in  this  study  is  also  from  the  1976 
U.S.  Standard  Atmosphere.  The  mixing  ratio  for  water-vapor  is  not 
constant  and  is  tabulated  in  the  reference  in  2  km  altitude  Increments. 
The  values  given  are  average  mid-latitude  concentrations  and  they 
Include  average  contributions  from  thunderstorms  and  other  stratospheric 
sources  above  16  km  (see  Ref  9:44  for  additional  detail).  The  water- 
vapor  number  densities  used  in  this  study  are  given  in  Appendix  I. 

In  the  infrared  wavelength  band  considered  in  this  study,  most  of 
the  atmospheric  scattering  of  electromagnetic  radiation  is  caused  by 
particulates  rather  than  by  gas  molecules.  Thus,  the  aerosol  model  used 
is  of  particular  Importance.  Three  different  models,  from  which  the 
aerosol  absorption  and  scattering  coefficients  can  be  calculated,  were 


examined  for  this  study.  The  model  by  Shettle  and  Fenn  (Ref  22)  was 
used  by  Wank  in  his  study.  An  earlier  (1964)  model  by  Elterman  (Ref  23) 
was  also  considered,  as  well  as  a  model  by  McClatchey  et  al.  published 
in  1972  (Ref  10). 

The  model  by  McClatchey  was  chosen  for  use  in  this  study.  Like  the 
Shettle  and  Fenn  model,  it  is  based  on  recent  experimental  measurements; 
but  it  has  the  added  benefit  that  it  presents  the  aerosol  extinction 
coefficients  in  tabular  form,  as  a  function  of  altitude,  and  at  selected 
wavelengths  that  span  the  visible  and  the  IR.  Coefficient  values  at 
other  wavelengths  can  be  linearly  interpolated  between  tabulated  values 
with  a  stated  error  of  less  than  t  8  percent.  The  interpolated  values 
at  2154.6050  cm“^  are  given  in  Appendix  C.  The  values  listed  are  for  a 
clear  atmosphere,  visibility  23  km  at  0.55  microns. 


Absorption  and  Scattering  in  the  Atmosphere 

When  electromagnetic  radiation  propagates  through  the  atmosphere, 
two  main  processes  are  in  action,  absorption  and  scattering.  These  two 
processes  must  be  considered  (1)  at  the  molecular  level,  where  the  wave¬ 
length  of  the  Incident  IR  radiation  is  much  larger  than  the  molecular 
diameter,  and  (2)  at  the  aerosol  level,  where  the  incident  wavelength 
and  particle  diameter  are  roughly  equivalent. 

Absorption  and  scattering  enter  into  the  Ildar  equation  at  two 
points.  First,  the  atmospheric  transmission  term  quantitatively  descri¬ 
bes  the  attenuation  of  the  propagating  signal  using  Beer's  law; 
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where  R  is  the  distance  propagated,  and  £  is  referred  to  as  the  extinc¬ 
tion  coefficient,  which  can  be  explicitly  written  as: 


c  -  e. 


«e4rT«a 


Mvi.ao'L/'*  ^  ^oLicoi.AA  (Sf) 

mnaiJ 


where  the  individual  contribution  of  each  process  has  been  recognized. 

The  two  aerosol  extinction  coefficients  are  usually  summed  and 
referred  to  as  the  Mie  extinction  coefficient,  . 

The  molecular  scattering  extinction  coefficient  is  generally 
referred  to  as  the  Rayleigh  scattering  extinction  coefficient  and 
notated  Cg  . 

The  molecular  absorption  extinction  coefficient  is  dependent  upon 
the  specific  absorption  effects  of  all  the  molecules  in  the  path,  and  is 
therefore  extremely  wavelength  sensitive.  In  this  study  the  line-by¬ 
line  contribution  of  five  molecular  species  are  considered ,  and  the 
extinction  coefficient  is  written  as  the  sum  of  five  terms  (see  equation 
31). 

The  second  term  in  the  lidar  equation  that  depends  on  the  absorp¬ 
tion  and  scattering  processes  is  the  volumetric  backscatterlng  coef¬ 
ficient,  jS  .  It  is  this  term  that  determines  how  much  incident  energy 
Is  scattered  back  toward  the  Ildar  receiver.  It  is  usually  written  as 
the  sum  of  two  terms : 


Is  due  to  the  Rayleigh,  or  molecular  scattering  process;  and 
Is  due  to  Mle  or  aerosol  scattering.  As  one  might  expect,  Alt  can  be 
related  to  ;  and  can  be  related  to  through  the  use  of  an 

angular  distribution  function. 

Rayleigh  Scattering  Extinction  Coefficient 

The  extinction  coefficient  due  to  Rayleigh  scattering  can  be  calcu¬ 
lated  from: 

where  Is  the  total  number  density  of  molecules ,  and  Is  the 

Rayleigh  scattering  cross-section. 

There  Is  a  quasi -empirical  equation  for  the  Rayleigh  scattering 
cross-section  that  yields  accurate  results  up  to  an  altitude  of  100  km 
(Ref  15:98); 

CA)  =  v.r<;  C A  /  o.st)  *  /o  (39) 

where  A  Is  the  wavelength  In  microns. 

Using  equation  ( 38)  and  the  total  atmospheric  molecular  number  den¬ 


sity  from  Appendix  B,  the  Rayleigh  scattering  extinction  coefficient  can 
be  calculated  as  a  function  of  altitude  for  a  wavelength  of  2154.6050 


The  volumetric  backscatter  coefficient  due  to  Rayleigh  scattering 


can  be  calculated  from  (Ref  2:22): 


^roT  (It) 

SO 

where  jT]  is  the  differential  scattering  cross-section  in  the  0  > 
180*  direction,  which  has  dimensions  cm^  sr“l,  and  N^qt  total 

number  density  of  all  molecules  in  cm'^. 

The  differential  scattering  cross-section  is  given  by  (Ref  2:23): 


s  _ 2  (  /  f  COS  e  J 

dSX.  TT 


(Vo) 


where  is  the  Rayleigh  scattering  cross-section  and  9  is  the  scatter 
angle  (  ®  *  0  is  in  the  forward  direction).  Therefore: 


d  r- 
dA 


3 

«  TT 


Then,  combining  equations  (37),  (39),  and  (41): 


(V/) 


(vc) 


where  the  factor  S/ST  is  recognized  as  the  angular  distribution  term 
noted  earlier. 


The  Rayleigh  volumetric  backscatterlng  coefficients  calculated  from 


equation  (42)  and  the  values  in  Appendix  D,  are  given  in  Appendix  E. 


Aerosol  Extinction  Coefficients 

Atmospheric  attenuation  due  to  absorption  and  scattering  by  par¬ 
ticulates  can  be  derived  theoretically  for  spherical  particles  from  a 
theory  developed  by  Mle  (1906).  'Hius  aerosol  scattering  is  referred  to 
as  Mle  scattering.  In  contrast  to  the  Rayleigh  scattering  phenomenon, 
Mie  scattering  is  a  great  deal  more  complex.  The  Mle  absorption  and 
scattering  extinction  coefficient  depends  on  the  number  density  and  size 
distribution  of  the  aerosols ,  as  well  as  on  their  complex  indices  of 
refraction.  Whereas,  the  Rayleigh  extinction  coefficient  could  be 
calculated  from  a  simple  empirical  expression,  the  Mie  extinction  coef¬ 
ficient  is  usually  tabulated  as  a  function  of  wavelength  and  altitude 
for  any  of  several  meteorological  conditions. 

The  Mie  extinction  coefficients  used  in  this  study  are  given  in 
Appendix  C.  As  stated  earlier,  is  the  sum  of  an  absorption  component 
and  a  scattering  component.  These  components  are  also  given  in  Appendix 
C. 


Aerosol  Volumetric  Backscattering  Cefficient 

The  aerosol  volumetric  backscattering  coefficient  is  given  by 
(Ref  2:16): 

-It  1 

■”  dJTL. )  ^ror  (ys) 

/eo 
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where  differential  backscaCterlng  cross-section,  this  time 

due  to  Hie  scattering;  and  is  the  total  particle  number  density. 

In  Che  section  on  the  Rayleigh  backscaCterlng  coefficient,  it  was 
shown  Chat  Ajf  and  £4  are  related  through  an  angular  distribution 
function.  In  analogous  fashion,  it  can  be  shown  Chat  (Ref  2:16): 

~  ^  AtAcsot,  CVV) 

SCATTtA 

where  f  is  the  distribution  function  of  interest  and  £  is 

«<4rra« 

the  scattering  component  of  .  Certainly  Che  temptation  here  is  Co 
assume  an  isotropic  distribution  and  set  fC»>B  t/i^ir  ^  since  the  actual 
distribution  function  is  not  analytically  known  for  particles  of 
arbitrary  size  and  surface  shape.  As  it  turns  out,  the  distribution 
function  is  not  isotropic,  but  experimental  measurements  indicate  Chat 
the  ratio  about  0.01-0.04  (Ref  15:96).  If  /(•)  is  asstimed  to 

be  i/lf  at  ■  180®,  a  sample  calculation  using  the  data  in  Appendix  C 
and  equation  (44)  yields  a  fin  ^ £n  ratio  of  0.04  at  sea-level. 
Therefore ,  using  an  angular  distribution  factor  of  /  /ttr  produces 
results  that  are  in  agreement  with  experimental  measurements. 

The  final  equation  for  fin  can  then  be  written  as : 

The  aerosol  volumetric  backscaCterlng  coefficients  calculated  from 
equation  (44)  and  Che  data  in  Appendix  C,  are  listed  in  Appendix  E. 

Also  listed  in  this  appendix  are  the  values  for  Che  total  volumetric 
backscaCterlng  coefficient  calculated  from  equation  (36). 


Molecular  Extinction  Coefficients 


The  atmospheric  extinction  coefficient  for  a  given  molecular  spe¬ 
cies  Is  given  by  (Ref  15:103): 


f.  =  /✓. 

t  t  * 


where  Is  the  number  density  of  the  absorbing  species ,  and  V,- 

Is  Its  absorption  cross-section.  The  absorption  cross-section  In  this 
equation  Is  dependent  upon  the  line-strength  of  the  of  the  particular 
absorption  line  of  Interest,  and  also  Its  line-width.  These  two  parame¬ 
ters  can  be  used  to  define  a  line-shape  function  which  gives  the  absorp¬ 
tion  cross-section  as  a  function  of  wavenumber. 

In  the  lower  atmosphere,  the  line-shape  function  Is  Lorentzlan  and 
takes  the  form  (Ref  10:12): 


C  v) 


.  /£]  _ _ 

V  (-K- n)*  H  * 


where 


Is  the  line  strength 

Is  the  line  half -width  at  half -maximum 

Is  the  wavenumber  of  the  Incident  radiation 

Is  the  wavenumber  of  the  absorption  line 
center 


The  Lorentzlan  line -shape  Is  a  product  of  collision  broadening,  and 
Is  valid  to  an  altitude  of  perhaps  20  km.  Very  high  In  the  atmosphere 
(above  30  km)  the  line-shape  Is  Gaussian,  which  Is  an  Indication  that 
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only  Doppler  broadening  effects  are  present.  In  the  region  between,  at 
20-30  inm  altitude,  the  line-shape  function  Is  a  convolution  of  the 
Lorentzlan  and  Gaussian  functions,  known  as  a  Volght  profile. 

This  study  assumes  that  the  absorption  line -shape  can  be 
represented,  throughout  the  altitudes  of  Interest,  by  a  simple 
Lorentzlan,  as  given  In  equation  (46).  The  cross-section  so  derived  Is 
approximately  correct  and  will  not  affect  the  accuracy  of  the  results, 
since  absorption  effects  above  20  km  will  be  small  In  any  case  (Ref 
10:12). 

The  line-strength  and  line-width  parameters  used  to  calculate  the 
absorption  cross-section  are  both  functions  of  altitude.  An  equation 
giving  the  temperature  dependence  of  line-strength  has  already  been 
developed  and  was  stated  earlier  as  equation  (26).  Using  this  equation 
and  the  temperature  data  in  Appendix  B,  the  line-strength  as  a  function 
of  altitude  can  be  calculated  from  its  sea-level  value. 

The  pressure  and  temperature  dependence  of  the  absorption  line- 
width  is  given  by  (Ref  10:12): 


<»<  <  T,  />  ) 


(.H7) 


where  ^  is  the  line  half-width  at  half-maximum  at 

sea-level 

Tg  is  the  sea-level  temperature 
P  /  is  the  atmospheric  pressure  ratio 


From  this  equation  and  the  data  in  Appendix  B,  the  absorption  line-width 
as  a  function  of  altitude  can  be  calculated  from  its  sea-level  value. 


Once  line -strength  and  line-width  are  known,  it  becomes  a  simple 
matter  of  using  equation  (46)  to  calculate  the  absorption  cross-section 
for  a  single  line  at  the  wavenumber  of  Interest.  This  approach  is  ade¬ 
quate  if  only  one  or  a  small  number  of  lines  contribute  to  the  total 
absorption  cross-section,  but  for  a  large  number  of  lines  the  number  of 
calculations  becomes  prohibitive. 

The  absorption  cross-section  at  altitude  can  be  approximated  from 
its  sea-level  value  if  the  wavenumber  at  which  the  cross-section  is  to 
be  calculated  is  not  in  the  central  peak  of  any  absorption  line.  If  this 
condition  holds,  Chen  Che  equation  for  absorption  cross-section  becomes: 

“  (i)  (V-V,)*  (V.) 

since  the  term  is  much  greater  than  for  wavenumbers  far 

from  y,  .  Thus,  if  the  sea-level  value  of  ^  is  given  by: 


(yf) 


and  its  value  at  altitude  by: 


and  finally,  using  equations  (26)  and  (47): 


= 


(rz) 


The  value  of  E"  used  In  the  approximation  should  be  associated  with  the 
absorption  line  that  contributes  most  to  the  total  cross -section. 

Using  the  equations  developed  above,  the  required  molecular  extinc¬ 
tion  coefficients  can  now  be  calculated. 


Extinction  Coefficients  for  CO.  The  molecular  extinction  coef¬ 
ficient  at  2154.6050  cm"^  due  to  nearby  CO  lines  (not  to  Include  the 
R(2)  line)  can  be  calculated  at  sea-level  using  equation  (45).  The 
total  absorption  cross-section  will  just  be  the  sum  of  the  Individual 
cross-sections  for  each  line  calculated  from  equation  (46). 

The  only  parameter  that  Is  needed  to  complete  these  calculations  is 
the  CO  line  half -width  at  half -maximum.  For  this  study,  the  experimen¬ 
tal  measurements  by  Bouanlch  and  Haeusler  will  be  used  (Ref  24:695). 

They  report  a  CO  absorption  line  half -width  of  about  0.07  cm“^  atm~^  for 
the  lower  J  numbered  transitions. 

The  calculated  sea-level  absorption  cross-section  values  for  the 
nearest  25  CO  lines  are  given  in  Appendix  F.  The  total  absorption 
cross-section  ,  is  calculated  to  be  2.595  x  10~21  cm2  at  sea- 

level. 

The  value  of  V",*-  ®  function  of  altitude  can  now  be  calcu¬ 

lated  from  equation  (52)  and  the  data  in  Appendix  B.  The  value  for  E” 
used  should  be  the  one  associated  with  the  R(3)  transition,  since  It 


contributes  most  to  the  total  sea-level  cross-section.  The  lover  state 
energy  for  this  transition  Is  ''alculated  from  the  Dunham  coefficients  as 
1104.65560  c  “1.  The  cross-section  values  as  a  function  of  altitude  are 
given  In  Appendix  G. 

Finally,  the  atmospheric  extinction  coefficient  ^eo- 
calculated  from  the  CO  number  density  data  In  Appendix  B,  and  the  cross- 
section  data  just  determined.  These  results  are  also  given  In  Appendix 

G. 

CO2  Absorption  Near  2154  cm~l.  Carbon  dioxide  is  a  linear  trlato- 
mlc  molecule  whose  vibration-rotation  bands  span  the  range  between 
approximately  500  cm“^  and  12,820  cm~l  (Ref  15:37).  The  CO2  molecule 
has  three  main  vibrational  frequencies:  ♦'i  ,  I'x  ,  and  Vj  .  Of  these, 
V)  Is  optically  Inactive  because  of  the  symmetry  of  its  dipole  moment 
during  vibration.  The  remaining  two  vibration  bands,  V,  and  Vj  ,  are 
centered  near  667  cm“^  and  2326  cm“l,  respectively.  The  main  Vj  band, 
along  with  the  combination  band  are  close  enough  to  the  primary 

Ildar  wavelength  to  require  further  examination. 

The  Vj  band  of  CO2  absorbs  very  heavily  In  the  wavelength  region 
2247  cm~l  to  2387  cm~l.  The  P(50)  line  for  this  band  is  located  at 
about  2242  cm~^,  and  the  line  structure  Is  weak  past  this  point  (Ref  25: 
5-95).  From  this  It  can  be  concluded  that  there  are  no  strong  lines 
from  this  band  present  at  the  primary  wavelength. 

The  y,  y'j,  combination  band  of  CO2  Is  centered  at  about  2079  cm~^ 
and  extends  from  about  2037  cm“^  to  2114  cm~^.  The  primary  Ildar  wave¬ 
length  at  2154  cm“^  lies  well  beyond  the  R(40)  line  of  this  band 


(Ref  25:5-95).  Therefore  it  can  be  concluded  that  any  absorption 
effects  due  to  the  CO2  molecule  at  the  primary  wavelength  considered  in 
this  study  will  be  small  and  can  be  neglected. 

N2O  Absorption  Near  2154  cm~^.  Nitrous  oxide  is  a  linear  asym¬ 
metric  molecule  with  three  main  vibration  frequencies:  V,  at  1285.6 
cm~^,  at  588.8  cm~^,  and  at  2223.5  cm~^.  Even  though  N2O 

has  many  combination  bands  and  hot  bands ,  they  do  not  contribute  much  to 
atmospheric  absorption. 

There  are  two  bands  of  N2O  near  the  primary  wavelength:  a  weak, 

(Vj  *  ♦i  )  -  ,  band  at  about  2212  cm“l  and  the  fundamental  Vj  band 

at  2223.5  cm~^.  The  Vj  band  has  strong  lines  extending  from  about 
2163  cm”l  to  2262  cm"l  (Ref  25:5-95),  The  weaker  band  extends  from 
about  2165  cm“l  to  2242  cm”l.  Therefore,  it  can  be  concluded  that  any 
absorption  effects  from  N2O  at  2154  cm“^  will  be  small  and  can  be 
neglected. 

Ozone  Absorption  Near  2154  cm~^^.  Ozone  is  an  asymmetric-top  mole¬ 
cule  with  three  main  vibration  frequencies.  All  three  are  active  in 
infrared  absorption.  Two  combination  bands  are  in  the  region  of 
Interest  to  this  study:  the  strong  band  centered  at  about  2105 

cm“l,  and  the  weaker  band  centered  at  about  2206  cm”^-. 

The  primary  wavelength  at  2154.6050  cm”^  is  situated  between  these 
two  bands.  The  stronger  of  the  two,  ^  Vj  ,  has  been  measured  at  high- 
resolution  (Ref  26:416).  The  extent  of  the  measured  lines  in  the  study 
ranged  from  2065.1  cm~^  to  2135.5  cm~^.  Here  again,  it  can  be  concluded 


that  any  absorption  effects  due  to  O3  will  be  small  and  can  be 


neglected 


Extinction  Coefficients  for  H2O.  The  wavelength  region  around  2154 
cm“^  contains  numerous  water-vapor  absorption  lines.  Water-vapor  has 
an  asymmetric-top  molecular  structure.  The  absorption  spectra  of  this 
molecule  are  extremely  complex  and  difficult  to  calculate.  Fortunately 
there  is  an  excellent  reference  available  by  Flaud  and  Camy-Peyret  (Ref 
27)  that  gives  the  line-positions,  line-intensities,  and  lower  state 
energy  levels  for  the  water-vapor  molecule  in  the  wavelength  region  of 
Interest  to  this  study.  These  values  were  calculated  at  a  pressure  of  1 
atm  and  a  temperature  of  296*K.  Appendix  H  lists  the  position  and  line- 
strength  for  each  of  58  water-vapor  lines  near  the  primary  wavelength. 
They  extend  from  about  2141  cm~l  to  2168  cm"^.  Only  absorption  lines 
from  the  H2^^0  Isotope  are  included,  since  it  comprises  0.99729  of 
atmospheric  water-vapor  (Ref  16:6). 

The  absorption  cross-section  at  2154.6050  cm“^  for  each  individual 
line  can  be  calculated  from  equation  (46).  The  value  for  the  H2O  line 
half -width  is  taken  to  be  0.032  cm~^  (Ref  16:20).  The  results  are  tabu¬ 
lated  in  Appendix  H. 

The  total  absorption  cross-section  of  H2O  at  sea-level  is  the  sum 
of  the  individual  cross-sections.  Its  value  is  calculated  to  be 
4.435  X  10-24  cm2  at  2154.6050  cm*!. 

The  total  absorption  cross-section  as  a  function  of  altitude  can  be 
calculated  from  this  sea-level  value  using  equation  (52).  The  H2O  line 
that  contributes  most  to  the  total  cross-section  is  the  one  at  2154.710 
cm”^.  The  lower  state  energy  for  this  line  is  given  as  1821.599  cm”^ 


(Ref  27:150).  Therefore,  this  value  should  be  used  for  E”  In  the  calcu¬ 


lations.  The  resulting  total  cross-sections  are  tabulated  in  Appendix  I 
The  atmospheric  extinction  coefficient  calcu¬ 

lated  from  equation  (45).  The  H2O  number  density  is  given  in  Appendix 
I,  as  is  the  total  absorption  cross-section.  The  resulting  extinction 
coefficients  are  also  listed  in  this  appendix. 


Calculation  of  ^ a 

The  atmospheric  transmission  term  ,  can  now  be  determined.  The 

value  of  7^  over  a  1  km  vertical  path  at  a  range  R  from  the  lidar 
system  is  given  by  equation  (33);  where  the  contributions  due  to  CO2, 
N2O,  and  O3  have  been  shown  to  be  negligible.  Equation  (33)  then 
becomes : 

(R)  ^  ^  ^  j  (53) 


where  the  tabulated  values  of  the  extinction  coefficients  are  available 

z 

in  the  relevant  appendlcles.  The  resulting  values  of  7^  (f/  at  a 
wavelength  of  2154.6050  cm~^  are  given  in  Appendix  J. 

The  value  of  7^  y  ,  that  is  the  two-way  transmission  over  a  path 
from  0  to  a  range  R,  is  then  given  by  equation  (32).  The  tabulated 
values  for  this  quantity  are  also  available  in  Appendix  J.  Thus,  one 

can  easily  determine  from  the  tables  that,  at  a  range  of  200  km 

a  < 

(sea-level),  the  value  ot  J 


is  0.6322 


z 

Calculation  of  Teo 

X 

The  transmission  term  7^*  ,  which  Is  due  solely  to  the  R(2)  line 

of  CO,  will  be  the  dominant  transmission  term  in  the  Ildar  equation.  As 
such,  its  value  must  be  known  precisely.  Since  only  a  single  absorption 
line  is  involved,  the  absorption  cross-section  at  the  primary  wavelength 
can  be  calculated  directly  at  each  altitude  using  equations  (26),  (46), 
and  (47).  The  line-strength  at  each  altitude  can  be  calculated  from  its 
sea-level  value  and  equation  (26),  where  E“  is  calculated  from  the 
Dunham  coefficients  as  1093.12109  cm~^.  Then  the  line-width  as  a  func¬ 
tion  of  altitude  can  be  determined  using  equation  (47).  Finally,  the 
cross-section  can  be  calculated  from  equation  (46). 

The  CO  R(2)  line  half-width  at  half -maximum  has  an  experimentally 
measured  value  of  0.07  cm”l  at  1  atm  (Ref  24:701). 

The  line-strengths,  line-widths,  and  absorption  cross-sections 
determined  from  the  above  equations  are  given  in  Appendix  K, 

Once  the  absorption  cross-section  is  known,  the  two-way  atmospheric 
transmission  through  a  1  km  vertical  path  at  a  range  R  can  be  calculated 
from  equation  (29).  The  results  are  given  in  Appendix  L. 

The  transmission  term  for  the  path  from  0  to  a  range  R  can 

then  be  determined  from  equation  (28).  These  results  are  also  given  in 
Appendix  L. 

The  Effects  of  Finite  Laser  Bandwidth 

The  lidar  equation  implicitly  assumes  that  the  transmitted  laser 
radiation  is  monochromatic.  This  assumption  is  widely  used  in  the 
literature,  and  is  approximately  true  for  most  low-power  laser  systems 


whose  line-widths  are  generally  100  MHz  (0.003  cm"^)  or  less.  Thus, 
their  spectral  width  is  much  less  than  the  width  of  a  typical  absorption 
line  at  sea-level,  which  is  about  0.10  cm~^.  However,  Korb  et  al.  has 
pointed  out  that  the  effects  of  finite  laser  bandwidth  can  cause  very 
large  errors  in  range -re solved  DIAL  systems,  and  therefore  they  must  be 
accurately  accounted  for  (Ref  28:78). 

This  study  assumes  the  use  of  high-power,  pulsed  CO2  TEA  laser 
sources  whose  nominal  linewldth  is  about  2-3  GHz  (0.07-0.10  cm"^). 
Therefore,  the  laser  linewldth  is  of  the  same  order  of  magnitude  as  the 
width  of  the  absorption  line  being  measured.  Clearly,  an  attempt  should 
be  made  to  Cake  Che  finite  bandwidth  of  the  laser  into  account. 

In  Che  Kbrb  study,  a  correction  was  applied  Co  experimentally 
measured  lidar  data  using  an  iterative  procedure  to  calculate  Che 
atmospheric  transmission  averaged  over  the  laser  line-width.  The 
average  transmission  was  given  as : 


^  are  the  wavenumber  limits  of  the  integration 

Z,  ,  *■  is  the  path  over  which  the  calculation  is 
performed 


Recall  equation  (3)  for  Che  received  power  at  Che  primary 


wavelength : 


t  I 

The  first  atmospheric  transmission  term,  J  ,  was  chosen  to  be  rela¬ 
tively  Insensitive  to  wavelength,  that  Is  It  does  not  contain  any 
contributions  from  the  central  peak  of  an  absorption  line.  Indeed,  the 
only  term  In  this  equation  that  Is  extremely  sensitive  to  wavelength  Is 
the  second  atmospheric  transmission  term,  7%,  ,  which  Is  just  due 

to  Che  single  R(2)  absorption  line  of  CO.  This  will  be  the  dominant 
transmission  term  since  most  of  the  absorption  will  be  due  to  this  one 
line. 

This  suggests  Chat  the  finite  bandwidth  of  the  laser  can  be  taken 

Into  account  by  considering  the  laser  Interaction  with  this  line.  That 

%  . 

Is,  instead  of  calculating  the  term  as  In  equation  (8),  use  a 

value  averaged  over  Che  spectral  width  of  the  laser  line.  Since  the 
parameters  of  the  absorption  line  are  functions  of  altitude,  the  average 
transmission  must  be  calculated  for  each  range-cell.  Then,  the  total 
transmission  for  the  path  of  interest  can  be  obtained  by  taking  a  pro¬ 
duct  of  Individual  transmissions,  as  before. 

Consider  the  average  transmission  function  for  one  range-cell: 

/■A* 


A  c-k)  e 
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£  Cv,  r  )  dr 


A  fv) 


In  order  to  accurately  calculate  the  average  transmission  for  the  cell, 
two  functions  must  be  known:  the  spectral  line  function  of  the  laser 
pulse,  kiv)  ;  and  the  shape  of  the  absorption  line,  Ccr) .  In  this 
study,  both  will  be  approximated  by  Lorentzian  functions.  The  absorp¬ 
tion  line  shape  Is  Lorentzian  for  altitudes  below  20  km,  which  is  where 
most  of  the  attenuation  takes  place,  so  the  approximation  should  be 
reasonable.  The  spectral  line  shape  of  a  CO2  TEA  laser  Is  also 
Lorentzian,  at  least  theoretically,  so  the  approximation  should  be  valid 
here,  also  (Ref  29:17). 

Therefore,  let: 


where  Is  the  extinction  coefficient  due  to  the  R(2)  line 

oc  Is  the  half -width  of  the,  absorption  line 

5  is  the  line-strength 


is  Che  CO  number  density 

Is  Che  center  wavenumber  for  Che  '*<(2)  line 


■y  is  Che  wavenumber  for  the  calculation 


Then,  Che  average  transmission  can  be  calculated  for  each  cell  by 
numerically  Integrating  the  following: 


-r.l  (K)  = 

The  path  Integral  has  been  replaced  by  a  simple  multiplication  because 
Che  extinction  coefficient  is  assumed  constant  across  Che  range  cell. 

The  average  transmission  values  calculated  from  equation  (58)  are 
tabulated  as  a  function  of  altitude  In  Appendix  M.  The  product  over  a 
total  path  0  to  range  R  Is  also  presented  in  Appendix  M. 

Equation  (58)  was  numerically  integrated  by  Simpson's  rule  on  a 
digital  computer.  The  absorption  line  parameters  used  were  as  given  in 
Appendix  K  for  the  CO  R(2)  line  centered  at  V,  «  2154.5960  cm“^.  The 
laser  line  was  assumed  to  be  centered  at  »i  -  2154.6050  cm“^,  with  a 
line  half-width  of  0.05  cm~^  which  corresponds  to  a  3  GHz  bandwidth 
laser  source.  The  integration  limits  were  from  2154.1  cm~^  to  2155.1 
cm“l.  These  limits  represent  the  points  where  the  laser  distribution 
function  decreases  to  about  IZ  of  Its  peak  value.  The  normalization 
constant  used  had  the  value  a  ■  0.339877562.  The  Integration  step-size 
was  selected  to  ensure  an  accuracy  of  four  significant  figures.  At  the 
higher  altitudes  where  the  absorption  line-width  Is  small,  the  step-size 
used  was  6.67  x  10“^  cm~l. 
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Figures  2  and  3  contrast  the  spectral  transmission  functions  at 
the  two  altitude  extremes,  sea-level  and  30  km,  respectively.  In  each 
figure  the  solid  line  is  the  spectral  distribution  function  for  the 
laser,  hcy)  ;  and  the  dashed  line  is  the  integrand  in  equation  (58), 
which  is  just  the  differential  transmission  3T / iY  .  The  dashed  line 
may  therefore  be  taken  as  an  indication  of  the  distribution  of 
transmitted  energy  after  it  has  twice  traversed  a  1  km  path  at  that 
altitude.  Note  that  at  an  altitude  of  30  km,  the  separation  between  the 
absorption  line  center  frequency  and  the  laser  center  frequency  is  noti¬ 
ceable  due  to  the  narrow  absorption  line -width. 

As  stated  earlier,  this  technique  is  not  usually  used  when 
addressing  problems  of  this  type.  Usually,  a  monochromatic  source  is 
assumed.  Therefore,  in  the  remainder  of  the  report,  the  results  will  be 
given  in  terms  of  (1)  a  monochromatic  source  assumption,  and  (2)  a 
finite  bandwidth  laser  source  assumption.  In  this  way  the  results 
obtained  using  the  two  methods  can  be  easily  compared. 

In  this  regard,  it  is  convenient  to  Introduce  at  this  point  an 
effective  cross-section  for  the  CO  R(2)  line  that  will  preserve  the 
integrity  of  the  equations.  Therefore,  let  b®  defined  such  that 

the  following  equation  is  true : 

(sf) 

The  effective  cross-sections  so  defined  are  tabulated  as  a  function  of 


altitude  in  Appendix  M. 
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Atnospheric  Transmission  Averaged 


AtitDspheric  Transmission  Averaged  Over  Laser  Line-Width  at  30  km. 


Minimum  Slgnal-to-Nolse  Requirements 


The  Ildar  equation  gives  the  return  signal  power  available  to  the 
signal  processing  portion  of  the  receiver.  The  ratio  of  the  return 
signal  power  level  to  the  system  noise  power  level  Is  termed  the  power 
signal -to-nolse  ratio  of  the  signal.  Detection  of  the  returning  signal 
from  a  given  range-cell  Implies  that  this  ratio  is  at  least  one. 
Therefore,  a  basic  requirement  for  system  performance  Is  that  both  the 
primary  and  the  reference  wavelength  slgnal-to-nolse  ratios  be  greater 
than  or  equal  to  unity. 

Another  consideration  for  a  DIAL  Ildar  system  Is  that  the  dif¬ 
ference  between  the  two  signals  must  be  measureable.  In  this  study, 
the  criterion  that  will  be  used  Is  to  require  that  the  difference  bet¬ 
ween  the  primary  and  reference  signals  be  greater  than  or  equal  to  the 
system  noise  level. 

These  slgnal-to-nolse  requirements  can  be  formally  stated  as: 


P,  2  P- 

F^-  F,  2  P»  <4.) 


(60) 

(.<•() 


where 


P,  Is  the  primary  return  signal  power 
P^  Is  the  reference  return  signal  power 
Is  the  receiver  system  noise  power 


These  requirements  can  be  written  only  In  terms  of  the  primary 
signal  power  by  considering  that  Pi  and  ?2  are  related  by  equation  (5). 

«  I 

Since  the  transmission  term  must  always  be  less  than  or  equal 

to  one,  equation  (61)  will  be  satisfied  If  equation  (60)  Is.  Therefore, 
using  the  definition  Si  »  Pi/Pn»  the  two  can  be  combined  to  require: 

^  (C3) 


Then,  after  solving  equation  (5)  for  P2  and  substituting  the 
result  Into  equation  (62),  equation  (62)  becomes: 


An  analysis  shows  that  the  right  side  of  this  Inequality  will  be  greater 

i  . 

than  or  equal  to  one  for  all  values  of  greater  than  or  equal 

to  V2*  Therefore,  the  two  signal-to-noise  requirements  can  finally  be 
restated  as : 


S,  t  i 
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(65) 

(d.A) 


These  two  Inequalities  determine  the  basic  performance  of  the  DIAL 


Ildar  system.  Satisfaction  of  the  first  condition  ensures  that  a 
measurement  can  be  made  over  long  paths  where  the  absorption  Is  high. 
Satisfaction  of  the  second  ensures  that  the  ratio  of  the  two  signals  Is 


measureable  over  paths  where  the  absorption  is  low  and  the  two  signals 
are  very  nearly  equal. 

Equations  (65)  and  (66)  are  formulated  only  in  terms  of  quantities 
evaluated  at  the  primary  Ildar  wavelength.  Therefore,  the  remainder  of 
the  study  need  only  consider  performance  at  this  wavelength.  If 
desired,  reference  wavelength  parameters  can  easily  be  obtained  using 
the  relationship  stated  In  equation  (5). 

The  signal -to-nolse  requirements  stated  above  can  be  reformulated 
Into  a  set  of  required  minimum  slgnal-to-nolse  ratios  as  a  function  of 
altitude  by  evaluating  the  term  each  altitude  (see  Appendices 

L  and  M).  The  results  are  shown  In  Figure  4  for  both  the  monochr jmatlc 
source  and  the  finite  bandwidth  source. 

The  results  indicate  that  very  large  primary  wavelength  slgnal-to- 
nolse  ratios  are  necessary  at  the  higher  altitudes  where  the  differen¬ 
tial  absorption  Is  small. 

An  example  here  might  serve  to  illuminate  the  situation.  Consider 

the  monochromatic  source  case  at  an  altitude  of  30  Icm  (R  «  170  km), 

1  , 

where  «  0.9972.  Suppose  that  the  system  noise  level  is 

-13 

10  Watts.  Then  from  Figure  4  the  minimum  required  slgnal-to-nolse 

ratio  Is  356.  This  equates  to  a  primary  return  signal  power  of 
-13 

356  X  10  Watts.  From  equation  (5)  the  reference  signal  power  is 

-13 

determined  to  be  357  x  10  Watts.  The  difference  between  the  primary 

-13 

and  reference  signals  Is  precisely  10  Watts,  which  Is  the  system 
noise  level.  Thus,  If  the  slgnal-to-nolse  ratio  was  anything  less  than 
356,  the  difference  between  the  two  received  signals  would  be  less  than 


the  system  noise  level,  and  therefore  the  two  would  be  Indistinguishable. 

This  result  Is  In  agreement  with  what  Is  generally  stated  for 
systems  using  the  DIAL  technique,  namely  that  with  current  technology  a 
signal  difference  (A  P/P)  of  the  order  of  1  percent  Is  detectable 


The  measurement  of  carbon  monoxide  in  the  atmosphere  using  the  dif¬ 
ferential  absorption  technique  requires  the  sensitive  detection  of  small 
amounts  of  infrared  laser  energy.  There  are  two  basic  techniques  that 
can  be  used:  direct  photon  detection,  and  heterodyne  detection.  Direct 
photon  detection  is  the  most  commonly  used  and  also  the  simpler  of  the 
two  to  implement.  Detection  is  achieved  by  directly  converting  incident 
photons  into  charge  carriers .  The  amount  of  current  produced  by  the 
device  is  then  measured  and  related  to  the  incident  photon  flux. 

Heterodyne  detection  is  generally  much  more  sensitive  than  direct 
detection,  but  requires  more  complicated  electronics  and  an  additional 
laser  source  for  the  down-conversion  process.  The  principle  involved  is 
identical  to  the  more  familiar  heterodyne  process  used  at  microwave  fre¬ 
quencies  ,  Its  implementation  in  the  optical  region  is  accomplished  by  com¬ 
bining  the  return  signal  with  a  laser  local  oscillator  signal  at  the 
surface  of  a  photomixer.  The  non-linear  mixing  process  converts  the 
incident  signal  to  a  lower  frequency  that  can  be  amplified  by  conven¬ 
tional  electronics. 

On  the  surface  it  appears  that  heterodyne  detection,  with  its 
superior  sensitivity,  should  be  chosen  over  direct  detection.  However, 
because  of  its  relative  simplicity,  direct  detection  would  be  preferred 
if  both  techniques  should  provide  adequate  detection  capability. 

Therefore,  this  study  will  address  both  direct  and  heterodyne  detection 
schemes . 


Direct  Detection 


Infrared  detectors  can  be  described  by  four  main  parameters : 
spectral  coverage,  response  time,  minimum  detectable  power,  and 
operating  temperature.  For  this  study,  the  detector  must  provide  ade¬ 
quate  sensitivity  at  4.64yu  and  have  a  response  time  of  less  than  6.67 
microseconds,  which  is  the  Integration  time  of  the  receiver.  The  detec¬ 
tor  which  most  nearly  meets  these  requirements  is  the  indium  antimonlde 
(InSb)  photovoltaic  infrared  detector  operating  at  a  temperature  of 

77®K.  Such  detectors  typically  have  a  normalized  detectivity  of  about 
11  1/9  -1 

10  cm  Hz  W  at  4.6/*  with  a  response  time  of  less  than  one  microse¬ 
cond  (Ref  25:11-78),  The  detector  parameters  used  In  this  study  are  given 
In  the  following  table.  They  are  for  a  commercially  available  Ford 
Aerospace  and  Communications  InSb  ISC-368  infrared  detector,  and  are 
identical  to  those  used  in  the  study  by  Wank  (Ref  4:60).  The  nor¬ 
malized  detectivity  quoted  in  the  table  is  for  an  angular  FOV  of  120 


Table  IV 

InSb  ISC-368  Infrared  Detector 


Sd 

(Sensitivity) 

2.0  A/W 

D* 

(Detectivity  120®  FOV) 

8.6  X  lOlO  cm  HzV2  W"! 

Ad 

(Detector  Area) 

2.827  X  10~3  cm2 

degrees .  The  detectivity  for  a  FOV  of  10"^  radians  can  be  calculated  by 
considering  a  quantity  called  D**,  which  is  related  to  D*,  but  is  inde¬ 
pendent  of  the  system  FOV.  For  a  flat-surfaced  detector  with  circular 


symmetry  (Ref  30:45): 


D**  =  s//^  (&7) 

where  6y^  is  the  half -angle  FOV  of  the  detector.  Using  this  equation 
and  the  values  in  the  table,  D**  *  7.448  x  10^0  cm  HzV2  sr"^  W“l.  Then 
D*  for  a  10“^  radian  FOV  will  be: 

0=0/ 


where  D**  is  given  above  and  6  is  5  x  10"^  radians.  Therefore,  the 
value  of  D*  corrected  for  FOV  is  1.49  x  10^^  cm  HzV2  W”^. 

The  power  slgnal-to-noise  ratio  for  a  photovoltaic  detector  and 
amplifier  combination  is  given  by  (Ref  31:69): 


S 

N 


where 


.  z 
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t'j  is  the  mean  signal  current 

is  the  mean  background  current 

is  the  mean  dark-current 

e  is  the  electron  charge 

S  is  the  electrical  bandwidth  of  the 

system 

7^  is  the  effective  temperature 

is  the  load  resistance 
/f  is  the  Boltzmann  constant 


The  first  three  terms  in  the  denominator  are  the  mean-squared  shot  noise 
current  from  the  signal,  background,  and  dark-current,  respectively. 

The  last  term  in  the  denominator  is  the  mean-squared  noise  current  due 


to  t::e  amplifier  electronics.  Depending  on  which  noise  source  gives  the 
largest  contribution,  the  detector  Is  said  to  be  photon,  background, 
dark -current ,  or  amplifier  noise  limited,  respectively. 

In  order  to  address  the  performance  of  the  detector,  the  background 
noise  power  must  first  be  known. 

Background  Noise  Power.  There  are  three  principle  noise  sources 
which  contribute  to  the  total  background  noise  that  will  be  seen  by  the 
detection  system.  They  are: 

1.  Radiance  from  the  earth 

2.  Radiance  from  the  sun 

3.  Radiance  from  the  sky 

The  thermal  radiance  from  the  earth  can  be  estimated  by  assuming 

that  the  earth  Is  a  300*  blackbody  with  an  emlsslvlty  of  0.85.  Then,  If 
the  receiving  telescope  Is  pointed  directly  toward  the  earth,  the 
background  power  will  be: 

e 

=  W(X)  AX  A^T;^  K  (7®) 


where 


WC^)  Is  the  spectral  radiance  of  the  source 

In  W  cm”^  sr"l  um"^ 

^  A  Is  the  wavelength  Interval  under 

consideration  In  microns 

Is  the  steradlan  FOV  of  the  system 

A  Is  the  area  of  the  receiving  optics 

In  cm^ 

7^  Is  the  one-way  atmospheric  transmission 

Is  the  efficiency  of  the  receiving 
optics 
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The  wavelength  Interval  is  dependent  upon  the  optical  filtering  used  in 
the  detection  system.  It  must  be  wide  enough  to  recover  all  of  the 
laser  energy,  in  this  case  about  3  GHz;  yet  small  enough  to  provide 
adequate  noise  rejection.  Commercially  available  Interference  filters 
can  be  used  here,  and  a  of  O.OI^m  is  reasonable  (Ref  25:7-127).  The 
wavelength  Interval  used  in  the  following  calculations  is  centered  at 
4.640//  and  extends  from  4.635^  to  4.64^  .  The  one-way  atmospheric 
transmission  can  be  obtained  from  Appendices  J,  L,  and  M.  For  a  monochro 
matlc  source  Tj^  *  3.417  x  lO"^,  For  a  finite  source  Tj^  *  1.783  x  10”^. 
Using  standard  blackbody  tables  (Ref  32),  W(A)AX  is  determined  to  be 
1.535  x  10~6  w  cm~2  sr“l.  Then  using  the  values  for  A^.,  K,  and 
from  Table  I: 

e  -/J 

P  Si  f.  704  X  to  w*rrs  (Monochromatic  Source)  (7») 

o 

s  S.O^S"*  to  wttrr%  (Finite  Bandwidth  Source)  (Ti> 

0 

The  thermal  radiance  from  the  sun  can  be  estimated  by  assuming  that 
the  sun  is  a  blackbody  source  at  5900*K  with  an  emisslvity  of  1.0.  The 
receiving  telescope  is  again  assumed  to  be  pointing  directly  at  the 
earth.  Energy  from  the  sun  must  traverse  the  atmosphere,  be  reflected 
from  the  surface  with  some  efficiency,  and  then  traverse  the  atmosphere 
once  more  to  enter  the  receiving  optics.  The  background  power  from  the 
sun  will  then  be : 
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A,  il 


^ov 


K  r-  7", 


^73) 


where 


Is  the  projected  area  of  the  sun 


A, 


s-t 


Is  the  distance  from  Che  sun  Co 
the  earth 


is  Che  reflectivity  of  Che  earth 

and  all  other  parameters  are  as  previously  defined.  The  factor  ir  '  comes 
from  Che  assumption  that  Che  earth  is  a  lambertlan  reflector.  Using  Che 
blackbody  cables  and  the  same  AX  as  before,  W(X)AX  ■■  4.302  x  10“^  W 
cm“2  sr”l.  Then,  using  the  values  from  Table  I;  along  with  an  As  =>  1.52 
X  10^®  m2  ,  Rg-E  *  1.49  X  10^1  m,  and  an  rg  =  0.15;  the  background  power 
from  Che  sun  is  estimated  to  be : 
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(Finite  Bandwidth  Source) 
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The  thermal  radiance  from  the  sky  can  be  estimated  from  experimen¬ 
tal  data  by  Bell  (Ref  33:1316).  He  estimates  that  for  a  clear  sky  at  4.6>i« 
and  at  an  angle  of  90®  (vertical),  W(X)  *  2.0  x  10”^  W  cm~2  sr”!  u“^^. 

Using  this  information,  the  background  power  from  the  sky  can  be  calcu¬ 
lated  from: 
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where  all  quantities  are  as  previously  defined.  Assuming  that  W(  )  is 
constant  over  the  small  wavelength  interval  dX ,  W(X)AX  ■  2.0  x  10“®  W 

cm~2  3r“l.  The  background  power  from  the  sky  is  then  estimated  to  be: 


The  total  background  power  can  now  be  estimated,  since: 
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The  results  are : 
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Noise  Equivalent  Power.  The  purpose  of  this  section  Is  to  calcu¬ 
late  the  noise  equivalent  power  (NEP)  for  the  direct  detection  case. 

The  performance  of  the  InSb  direct  photon  detector  Is  dependent  upon  the 
contributions  of  the  various  noise  terms  In  the  denominator  of  equation 
(69).  One  noise  source,  that  associated  with  the  background  viewed  by 
the  detector,  has  already  been  estimated.  It  remains,  then,  to  estimate 
the  noise  contributions  due  to  signal  current,  dark -current ,  and  the 
amplifier  electronics. 

The  detector  dark-current  is  related  to  the  detector  normalized 
detectivity  by  the  following  (Ref  34:1648): 


(81) 


Therefore,  the  dark-current  will  be: 

“  2  «  D**- 

The  dark-current  calculated  from  this  equation  and  the  data  In  Table  IV 
Is  1.59  X  10~1^  A. 

The  mean  current  associated  with  the  background  signal  can  be 
calculated  from  the  product  of  the  background  noise  power  and  the  detec¬ 
tor  sensitivity.  The  mean  background  current  Is  about  8.0  x  10”^^  A. 
Since  the  background  noise  contribution  Is  three  orders  of  magnitude 
larger  than  that  due  to  dark -current ,  the  dark-current  noise  can  be 
neglected  In  the  following  calculations. 

The  mean-squared  noise  current  due  to  the  amplifier  and  Its  asso¬ 
ciated  electronics  Is  given  by  the  term  4kTeB/RL,  where:  Tg  is  the 
effective  temperature  of  the  detector/amplifier  combination,  Is  the 
load  resistance,  and  B  Is  the  electrical  bandwidth  of  the  system.  As 
stated  earlier,  this  bandwidth  must  be  adequate  to  support  the  6.67 
microsecond  integration  time  of  the  receiver.  For  this  study,  a  1  MHz 
bandwidth  is  selected.  This  should  provide  the  system  with  a  conser¬ 
vatively  chosen  response  time  of  about  1  microsecond,  as  well  as  ade¬ 
quate  noise  rejection. 

A  reasonable  value  for  the  effective  temperature  of  a  commercially 
available  low-noise  amplifier  chain  is  Tg  ■  250*K  (Ref  15:323). 

The  value  of  Rl  Is  dependent  upon  the  circuit  capacitance  and  the 
electrical  bandwidth  of  the  system.  A  typical  value  for  the  circuit 
capacitance  of  a  photodiode  detector  coupled  with  an  amplifier  Is 
Cj  ■  10  picofarads.  Then,  since  B  *  1/2  RlCj,  the  load  resistor  must 


have  a  value  of  about  10  k a  in  order  to  Insure  a  bandwidth  capability 
of  I  MHz. 

Substitution  of  these  values  into  the  equation  for  the  amplifier 
mean-squared  noise  current  yields  a  value  of  1.380  x  10“^®  a2. 

Equation  (69)  can  now  be  written: 
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Then,  incorporating  the  detector  sensitivity  into  the  equation  with  Ig 
Ps  Sd  and  ifl  -  Pg  Sq  : 
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where  Pg  is  the  signal  power,  and  Pg  is  the  background  noise  power. 

The  system  noise  equivalent  power  is  that  value  of  signal  power 
which  will  produce  a  power  signal-to-noise  ratio  of  unity.  Equation 
(84)  can  then  be  solved  to  yield: 
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Since  all  of  the  required  parameters  are  known  or  have  been  esti¬ 
mated,  the  system  NEP  for  both  the  monochromatic  laser  source  and  the 
finite  bandwidth  laser  source  can  be  calculated.  The  results  are: 
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(Monochromatic  Source) 


=  S’,  87r  X  /o'"'*  (Finite  Bandwidth  Source)  (87) 
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The  results  are  the  same  for  both  cases  because  the  dominant  term 
in  equation  (85)  Is  the  one  associated  with  the  amplifier  noise. 
Therefore,  In  the  direct  detection  case,  the  system  Is  shown  to  be 
amplifier  noise  limited. 

Required  Laser  Power.  Once  the  noise  equivalent  power  of  the 
system  Is  known,  the  minimum  slgnal-to -noise  requirements  developed 
earlier  can  be  addressed  to  yield  a  performance  assessment  of  the  InSb 
detector.  Consider  that  for  any  given  range  cell: 


=  (5/v)^ 


where  (S/N)min  is  the  required  minimum  slgnal-to'nolse ,  and  Pnep  is 
the  system  NEP.  Then  the  basic  Ildar  equation  can  be  solved  for  the  peak 
laser  power  required  to  achieve  that  minimum  signal -to-nolse  ratio.  That 
Is,  let: 

- - - - -  (8^) 

*  ^  (fZ-;  ^  r;/,  r„7^ 

The  required  laser  power  at  each  altitude  Is  shown  In  Figure  5. 
Results  are  displayed  for  both  the  monochromatic  source  and  the  finite 
bandwidth  source.  The  required  laser  energy  (assuming  a  100  ns  pulse) 
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is  also  shown 
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Heterodyne  Detection 

Heterodyne  detection  Is  the  most  sensitive  technique  for  detecting 
low-power  laser  radiation  in  the  infrared  portion  of  the  spectrum. 
However,  the  Increase  in  sensitivity  is  not  without  its  drawbacks,  since 
heterodyne  detection  imposes  stringent  requirements  on  the  system. 

These  requirements  are  basically  in  the  form  of  spatial  and  temporal 
coherence  considerations  that  must  be  met  in  order  to  achieve  the  full 
benefits  of  the  heterodyne  process. 

The  basic  heterodyne  detection  scheme  is  as  follows  (see  Figure  6). 
Laser  radiation  is  transmitted  into  the  atmosphere  and  scattered  off  of 
molecules  and  particulates  in  a  particular  range  cell.  The  returning 
photon  flux  is  collected  by  the  receiving  optics  and  applied  to  the  sur¬ 
face  of  a  photomixer,  where  it  is  combined  with  the  signal  from  a 
coherent  laser  local  oscillator  (L.O. ).  The  local  oscillator  frequency 
is  slightly  offset  from  the  signal  frequency,  and  their  difference  is 
referred  to  as  the  Intermediate  frequency  (IF)  of  the  process.  The  IF 
signal  contains  the  same  information  as  the  incoming  optical  signal,  but 
at  a  frequency  that  is  easily  handled  by  conventional  electronics.  The 
IF  or  difference  frequency  is  a  product  of  the  non-linear  photomixing 
process,  and  the  efficiency  of  this  process  determines,  in  large  part, 
the  sensitivity  of  the  detector. 

When  the  signal  wavefront  and  the  local  oscillator  wavefront  arrive 


at  the  surface  of  the  photomlxer,  they  must  be  spatially  coherent  in 
order  to  achieve  aximum  conversion  efficiency.  This  spatial  coherence 


requirement  can  be  stated  as  (Ref  35:37-47): 
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where  fi  Is  the  radian  angle  between  the  signal 

and  local  oscillator  wavefronts 

A  Is  the  wavelength  of  the  Incoming 
signal 

Jl_  Is  the  distance  over  which  the  photo- 

mlxlng  takes  place,  l.e.  the  diameter 
of  the  detector 


This  requirement  may  be  Interpreted  as  a  basic  limitation  on  the  system 
FOV,  where  —  iP  If  the  L.O.  signal  Is  aligned  with  the  optical 

axis  of  the  system.  For  a  wavelength  of  4,64^  and  a  detector  diameter 
of  1.0  X  10~3  m  (Ref  35:37-47),  fi  must  be  less  than  1.48  x  lO"^ 
radians.  This  leads  to  a  system  FOV  requirement  of  less  than  2.95  x  10"^ 
radians,  which  Is  easily  met  by  the  system  under  consideration. 

There  Is  one  other  spatial  coherence  requirement  that  must  be 
addressed.  The  signal  from  a  given  scatter  volume  Is  constantly  changing 
as  the  molecules  and  particles  move  about.  It  can  be  shown  that  the 
system  receiving  aperture  must  not  be  larger  than  the  mean  grain  size  In 
the  scatter  volume.  For  the  system  under  consideration,  the  requirement 
Is  that  the  area  of  the  receiving  optics  be  less  than  or  equal  to  the 
area  of  the  transmit  optics  (Ref  2:67).  Since  the  system  in  this  study 
uses  the  same  telescope  for  transmitting  and  receiving,  this  requirement 
Is  also  met. 

There  are  also  temporal  coherence  criteria  that  must  be  met.  They 


are  addressed  later  in  the  study.  At  this  point.  It  Is  convenient  to 


state  that  all  the  basic  requirements  for  heterodyne  detection  can  be 
satisfied  by  the  system  under  consideration  in  this  study. 


Noise  Equivalent  Power.  The  basic  limitation  on  the  sensitivity  of  the 
direct  photon  detector  was  determined  to  be  the  amplifier  noise.  In 
heterodyne  detection,  amplifier  noise  and  the  other  usual  noise  sources 
are  not  a  problem  because  of  the  manner  in  which  the  signal  is  pro¬ 
cessed.  For  heterodyne  detection,  the  power  slgnal-to-noise  ratio  can 
be  shown  to  be  (Ref  35:37-37): 
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that  is,  the  local  oscillator  signal  raises  both  the  signal  level  and 
the  noise  level.  Therefore,  amplifier  noise  and  other  noise  sources  can 
be  overcome  simply  by  raising  the  L.O.  signal  level. 

The  mean-squared  shot-noise  current  due  to  the  L.O.  signal  is: 
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where  is  the  quantum  efficiency  of  the  photomixer.  The  quantity  in 
brackets  is  just  the  mean  L.O.  signal  current.  If  the  L.O,  is  to  be  the 
dominant  noise  source,  then: 
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The  shot-noise  current  due  to  signal  and  background  were  determined  to 
be  much  less  than  that  due  to  the  amplifier  noise  in  the  section  on 
direct  detection.  Therefore,  the  first  term  can  be  neglected,  and: 
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is  all  that  is  required,  la  practical  terms,  the  required  L.O.  power  is 
usually  about  1  or  2  mW  (Ref  15:314).  Therefore,  the  requirement  that  the 
shot-noise  from  the  L.O.  be  the  dominant  noise  term  is  easily  attained. 

Assuming  that  the  local  oscillator  power  is  sufficient,  the  slgnal- 
to-noise  ratio  can  now  be  written: 
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where  is  the  mean-squared  signal  current,  which  can  be  obtained 

from: 
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Then,  the  slgnal-to-noise  ratio  becomes: 
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(97) 


Equation  (97)  is  the  form  in  which  the  heterodyne  signal-to-noise  ratio 
is  usually  written.  From  this  equation  it  is  easy  to  see  that  the 
heterodyne  noise  equivalent  power  is  given  by; 


7.  h  -v  B 


The  electrical  bandwidth  used  In  the  heterodyne  detection  process 
is  generally  matched  to  the  reciprocal  of  the  laser  output  pulse  dura¬ 
tion  (Ref  28:83).  For  a  100  nanosecond  pulse,  the  bandwidth  is  10^  Hz. 
Using  this  figure  for  B,  and  a  quantum  efficiency  of  0.5,  the  heterodyne 
NEP  is  calculated  to  be : 
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Heterodyne  Signal-to-Noise  Considerations.  The  slgnal-to-nolse 
ratio  given  in  equation  (97)  represents  near  quantum  limited  detection 
of  the  incoming  signal.  The  large  slgnal-to-nolse  ratios  Implied  by 
this  expression  are  attainable  only  if  the  received  photon  flux  is 
totally  coherent  during  the  receiver  Integration  time.  In  practice,  the 
return  signal  will  be  coherent  only  over  a  time  period  equivalent  to  the 
laser  pulse  length.  The  system  in  this  study  is  designed  to  integrate 
the  return  signal  for  much  longer  than  this  in  order  to  achieve  adequate 
sensitivity  over  very  long  ranges.  This  temporal  coherence  problem  is 
referred  to  in  the  literature  as  "speckle  noise”,  and  is  usually  handled 
by  considering  the  statistical  nature  of  the  return  signal.  The  deriva¬ 
tion  used  in  this  study  follows  closely  that  recommended  by  ^e  (Ref  36: 
3862). 

Consider  a  system  with  an  integration  time  that  is  long  compared  to 


the  pulse  duration.  Then,  the  return  signal  may  be  thought  of  as  being 
composed  of  a  series  of  partially  coherent  pulses  from  various  scatter 
points  in  the  range  cell  being  measured.  If  each  such  pulse  is  considered 


This  shows  that  for  any  given  sample,  the  achievable  slgnal-to-nolse 
ratio  Is  less  than  one,  so  long  as  the  photon  count  Is  sufficiently 


great.  Note  that  the  slgnal-to-nolse  ratio  Is  proportional  to  the 
received  photon  flux  only  In  the  small  signal  regime.  This  leads  to  a 
concept  of  saturated  slgnal-to-nolse  In  the  heterod5me  process  which  can 
only  be  Improved  by  Increasing  the  number  of  samples  m  ,  or  by 
averaging  a  number  of  laser  shots.  Then,  the  final  heterodyne  slgnal- 
to-nolse  ratio  can  be  written: 

(1)^  = 

or : 
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where  m  Is  the  number  of  samples,  and  M  Is  the  number  of  laser  shots 
averaged . 

The  requirement  that  »  I  Is  actually  quite  conservative. 

At  a  wavelength  of  jji  one  photon  contains  4.3  x  10“20  joules,  so 
over  one  pulse  period  the  expected  photon  count  at  the  system  NEP  Is 
about  four.  Therefore,  the  system  NEP  remains  a  valid  Indicator  of  sen¬ 
sitivity. 

Required  Laser  Power.  In  the  heterodyne  detection  case,  the 
attainable  slgnal-to-nolse  ratio  Is  Independent  of  the  laser  transmitted 
power,  so  long  as  a  sufficient  photon  count  Is  maintained.  As  shown  In 
the  last  section,  the  system  NEP  Is  an  accurate  Indicator  of  photon 
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count.  Therefore,  Che  single-shot  laser  power  required  to  meet  the 
signal-Co-noise  criteria  developed  earlier  can  be  calculated  from: 
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or : 


This  equation  is  valid  until  the  minimum  signal-to-noise  required 
exceeds  chat  achievable  with  heterodyne  detection  of  a  single  laser 
shot.  Past  this  point,  multiple  laser  shots  must  be  averaged  to  make 
the  measurement. 

For  the  system  under  consideration,  the  value  of  (S/N)n  is  given 
by  equation  (104),  where  M=-l  and  m=»66.7  .  The  result  is  (S/N)h  *  5.77 
Then,  from  Figure  4,  Che  system  can  make  single-shot  measurements  from 
sea-level  to  an  altitude  of  20  km  for  the  monochromatic  source  assump¬ 
tion,  and  to  an  altitude  of  17  km  for  the  finite  bandwidth  source. 

The  required  laser  transmitter  power  levels  for  these  altitudes  are 
shown  in  Figure  7. 


Error  Analysis 

Thus  far,  the  criteria  developed  to  assess  the  performance  of  Che 
lidar  system  have  been  designed  to  insure  that  Che  primary  and  reference 
signals,  and  their  difference,  will  be  measureable.  But  it  is  not 
enough  just  Co  be  able  Co  make  the  measurement.  The  purpose  of  this 


^  WAr-rs  ) 


Alt/  tups.  (km) 


Figure  7.  Laser  Power  Required  -  Heterodyne  Detection 
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section  Is  to  determine  the  accuracy  with  which  the  carbon  monoxide 
number  density  can  be  measured. 

Recall  that  the  equation  for  the  measured  number  density  is: 
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The  uncertainty  In  the  measurement  can  be  obtained,  assuming  Independent 
error  contributions,  from  the  standard  equation  for  propagation  of 
errors ; 
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Then,  assuming  that  the  average  absorption  cross-section  across  the 
range  cell  Is  known  exactly,  the  fractional  uncertainty  in  the  measure¬ 
ment  can  be  shown  to  be : 


A  lower  bound  on  the  expected  measurement  error  can  be  determined 
by  assuming  that  the  uncertainty  In  each  Individual  power  measurement  Is 
equal  to  the  system  noise  level.  Also,  the  equation's  dependence  on  the 


off-line  measurements  can  be  eliminated  by  using  the  relation  P2(K)  * 


Pl(R)  /  T,l  . 

Then  the  equation  becomes : 
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where 

•S/^a 

are  the  primary  wavelength 
signal-to-noise  ratios  at 
ranges  and  R2,  respectively 

is  the  two-way  atmospheric 
transmission  to  range  R^ 

is  the  two-way  atmospheric 
transmission  to  range  R2 

The  above  equation  is  exact,  but  can  become  tedious  to  work  with 
because  of  the  logarithmic  dependence  of  the  atmospheric  transmission. 
The  equation  can  be  further  simplified  by  using  a  few  approximations. 
Consider  that  the  two  signal-to-noise  ratios  in  the  equation  will  be 
approximately  equal,  since  they  are  for  adjacent  range  cells.  Also,  the 
two  transmission  terms  each  must  have  a  value  between  0  and  1.  Now,  if 
Sj  ■  $2  *  S  and  Im  ■  I.  “  0,  the  square-root  term  goes  toJz/S  . 
But,  If  Si  •  $2  ^  S  and  TrlL  ~r,li  -  1,  then  the  square-root  term 

goes  to  z/s  .  Finally,  consider  that  the  ratio  T„  /  is  just 

»»  », 

(X«).  The  above  considerations  suggest  that,  to  a  good  approxima¬ 
tion,  the  fractional  uncertainty  in  the  measurement  at  a  range  R  can  be 


written: 


where 


s 


is  the  power  signal-to-noise  at 
range  R 

is  the  two-way  atmospheric 
transmission  through  a  1  km 
vertical  path  at  range  R 


This  result  is  intuitively  satisfying,  since  it  takes  into  account 
the  two  limiting  factors  in  any  such  measurement.  First,  the  signal-to- 
noise  ratio  must  be  sufficiently  high;  and  second,  the  absorption  in  the 
range  cell  must  be  great  enough  that  the  difference  in  transmission  bet¬ 
ween  adjacent  cells  can  be  accurately  determined. 


VII.  Performance  Analysis  /  I  Joule  Laser  with  Heterodyne  Detection 


The  final  system  configuration  selected  for  a  performance  analysis 
utilizes  a  I  joule/pulse  CO2  TEA  laser  at  each  wavelength,  and  the 
heterodyne  detection  technique.  A  Ildar  system  with  these  charac¬ 
teristics  Is  well  within  the  state-of-the-art  and  most  of  the  components 
are  commercially  available.  A  summary  of  the  major  system  parameters  Is 
given  In  the  following  table. 


Table  V 

Lldar  System  Parameters 


Primary  (on-line)  wavelength 

2154.6050  cm-1 

Reference  wavelength 

2143.7674  cm'l 

Spectral  llnewldth 

0. 10  cm"^ 

Laser  repetition  rate 

100  Hz 

laser  pulse  duration 

100  ns 

System  NEP 

1.714  X  10-12  w 

Operational  altitude 

200  km 

Assuming  a  100  Hz  repetition  rate  and  a  20%  overall  efficiency  for 
both  the  primary  and  reference  lasers,  the  total  average  power  required 
for  the  system  should  be  on  the  order  of  I  kilowatt.  This  is  well 
within  the  limits  set  for  the  space  shuttle  (Ref  2:88). 
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The  return  signal  power  can  be  calculated  using  the  basic  Ildar 
equation  and  the  atmospheric  data  found  In  the  various  appendices.  The 
value  for  Pq  can  be  calculated  by  dividing  the  laser  pulse  energy  by  the 

pulse  duration.  For  this  system  the  peak  pulse  power  Is  10^  Watts.  The 

*  . 

atmospheric  transmission  term  ~^eo  I  should  be  obtained  from  Appendix  M, 

It 

since  the  laser  source  Is  assumed  to  have  an  output  llnewldth  of  about  3 
GHz . 

The  results  of  these  calculations  are  shown  as  a  function  of  alti¬ 
tude  above  sea-level  In  Figure  8.  As  can  be  seen  from  the  figure,  the 
return  signal  power  levels  are  well  above  the  system  NEP  across  the 
entire  altitude  spectrum. 


Even  though  the  return  signal  power  levels  are  well  above  noise, 
the  single-shot  slgnal-to-nolse  ratio  for  the  system  Is  limited  to  5.77 
by  the  nature  of  the  heterodyne  process.  This  limits  the  altitude 
regime  where  measurements  can  be  made  to  from  sea-level  to  about  13  km. 
In  order  to  meet  the  larger  slgnal-to-nolse  requirements  needed  at  the 
higher  altitudes,  successive  laser  shots  mu  it  be  averaged.  The  Improve¬ 
ment  In  slgnal-to-nolse  goes  as  the  square-root  of  the  number  of  pulses 
averaged,  as  shown  In  equation  (104).  However,  shot-averaging  will  also 
effect  the  spatial  resolution  of  the  measurements.  At  a  200  km  orbital 
altitude  the  velocity  of  the  space  shuttle  Is  about  7  km/s  (Ref  2:78). 


Thus,  if  100  pulses  are  averaged  at  a  laser  repetition  rate  of  100  Hz, 


the  horizontal  resolution  is  Increased  from  about  70  meters  to  7  km. 
Considering  the  large  volume  of  the  earth's  atmosphere  and  the  vertical 
resolution  of  1  km,  a  horizontal  resolution  of  7  km  Is  easily  accom¬ 
modated. 

Assuming  that  100  laser  shots  are  averaged  for  each  measurement, 
the  heterodyne  signal-to-noise  ratio  is  increased  by  a  factor  of  ten  to 
57.7  .  Then  from  Figure  4  the  system  should  be  able  to  make  measure¬ 
ments  from  sea-level  to  an  altitude  of  about  27  km. 

The  measurement  accuracy  is  given  by  equation  (110),  and  the  results 
are  shown  as  a  function  of  altitude  in  Figure  9.  They  Indicate  that, 
although  the  system  can  make  measurements  to  27  km,  the  accuracy 
will  not  be  good  enough  to  provide  useful  data  above  an  altitude  of 
perhaps  15  km  (50%  error  point).  The  reader  is  cautioned  that  the 
errors  shown  are  approximate  and  should  only  be  used  as  a  general  indi¬ 
cation  of  system  performance. 

Performance  Envelope 

Thus  far,  the  system  performance  has  been  measured  against  the  a 
priori  vertical  profile  of  CO  concentration  established  at  the  outset  of 
the  study.  The  atmospheric  transmission  parameters  were  all  calculated 
based  on  this  particular  profile.  Suppose  the  question  is  posed:  how 
much  or  how  little  CO  can  be  detected  at  any  particular  altitude?  This 
question  can  be  answered,  but  only  within  certain  constraints.  First  of 
all,  the  measured  concentration  at  a  particular  altitude  depends  on  the 


CO  concentration  in  the  intervening  layers.  That  is,  the  system  obtains 


a  measurement  by  comparing  the  returns  from  two  adjacent  range  cells; 
and  these  returns  are  a  function  of  the  CO  density  along  the  entire  ver¬ 
tical  path  from  the  shuttle  to  the  range  cell  of  interest.  The  approach 
that  will  be  used  here,  is  that  the  CO  mixing  ratio  is  fixed  at  0.075 
ppmv  at  all  altitudes  except  in  the  range  cell  being  analyzed.  This 
corresponds  to  an  assumption  that  a  large  gradient  in  the  CO  density 
profile  exists  at  a  certain  altitude.  The  system  capability  to  measure 
such  a  gradient  is  then  of  interest. 

The  maximum  amount  of  CO  that  can  be  measured  will  be  determined  by 
"pump  depletion";  that  is,  the  CO  number  density  is  so  great  that  most 
of  the  laser  power  is  absorbed,  and  the  return  signal  is  reduced  to  the 
system  noise  level.  Then  the  maximum  amount  of  CO  that  can  be  detected 
can  be  calculated  from  the  following: 


P,  (K)  = 

or : 


and  then,  since: 


NtP 


cm) 


(ift) 


I. 


A-l 


Tco  CK) 


(Hi) 


one  finally  obtains: 


(^‘)  (A.)  7-;/^ 
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Once  the  transmission  in  the  range  cell  is  calculated,  it  is  a  simple 


matter  to  obtain  the  maximum  number  density: 

HA*  1  r  X  1 

A/  -  -  - — —  Jn  / 

X  rCR'i  aK  L  J 

where  the  argument  of  the  natural  logarithm  is  given  by  equation  (114) 
above .  Here ,  AK  -  1  km  and  is  the  effective  cross-section  given 

in  Appendix  M.  An  effective  cross-section  must  be  used  because  of  the 
finite  bandwidth  of  the  laser.  It  should  be  noted  that  this  analysis 

X 

neglects  the  small  change  in  due  to  the  additional  CO. 

The  minimum  value  of  CO  that  can  be  detected  depends  on  the  dif¬ 
ference  in  transmission  between  two  adjacent  range  cells  that  can  be 
accurately  measured.  Recall  that  the  fractional  uncertainty  in  the 
number  density  measurement  is  given  approximately  by: 


^  ^  _ i _ 

then,  for  any  given  uncertainty,  the  transmission  across  the  range  cell 
must  be : 

T,]  iR)  -  exP  < - ; — C//4 

A  little  thought  will  show  that  this  value  of  Te.CO  represents  the 
maximum  transmission  (hence,  the  minimum  absorption)  that  can  be 
measured  with  the  stated  uncertainty.  Then,  the  minimum  number  density 


that  can  be  measured  in  the  range  cell  becomes : 


CO 


Z  AK  T-f*) 


L  'CO  \  f  J 


where  the  argument  of  the  natural  logarithm  Is  given  by  equation  (116) 
above.  This  leads  to: 


N 


CO 


_ i _ 

2  A)?  s  CdA/,, 


(n8) 


The  minimum  number  density  that  can  be  detected  Is  then  a  function  of 
how  accurate  the  measurements  must  be. 

Using  equations  (115)  and  (118),  a  performance  envelope  can  be 
constructed  for  the  system  under  consideration.  Figure  10  shows  the 
expected  performance  as  a  function  of  altitude.  The  standard 
atmosphere  profile  of  CO  number  density  Is  shown  for  reference.  Once 
again,  note  that  the  minimum  and  maximum  number  densities  that  can  be 
measured  are  based  on  the  assumption  that  a  standard  CO  density  profile 
exists  In  the  atmosphere  between  the  shuttle  and  the  range  cell  of 
Interest. 


Detection  of  "Hot"  CO 

Although  the  primary  function  of  a  Ildar  carbon  monoxide  detection 
system  In  space  would  be  to  measure  and  characterize  the  atmospheric  CO 
budget ,  a  secondary  function  might  be  the  detection  of  carbon  monoxide 
"hot”  spots  on  or  near  the  surface  of  the  earth.  Such  "hot"  CO 


emissions  might  be  from  a  factory  stack  or  perhaps  from  an  area  where 


the  automobile  density  Is  extremely  high.  Another  source  might  be  the 
exhaust  plumes  from  gas-turbine  or  rocket  engines. 

These  emissions  would  differ  from  the  standard  atmospheric  CO  con¬ 
tent  in  that  the  temperature  of  the  molecules  would  be  higher  than 
ambient.  The  range  of  temperatures  considered  in  this  study  is  from 
400*K  to  1000*K.  This  range  should  be  adequate  to  describe  the  pheno¬ 
mena  of  Interest. 

A  detailed  analysis  of  system  capabilities  in  this  area  is  beyond 
the  scope  of  this  thesis.  However,  a  general  idea  of  system  performance 
can  be  obtained  if  several  assumptions  are  made : 

1.  That  the  laser  illuminates  a  uniform  volume  of  "hot"  CO  . 

2.  That  the  “hot”  CO  has  diffused  to  the  point  where  the 
ambient  pressure  is  valid. 

3.  That  the  "hot”  CO  is  confined  to  an  altitude  of  less  than 
2  km. 

A,  That  the  "hot"  CO  is  present  only  in  the  range  cell  under 
consideration. 

If  these  assumptions  are  made,  then  the  maximum  "hot”  CO  number 
density  that  can  be  detected  will  be  that  value  of  which  will 

reduce  the  return  signal  power  to  the  system  noise  level.  The  return 
signal  without  the  "hot"  CO  has  already  been  calculated  and  is 
available  from  Figure  8.  Then  the  return  signal  power  with  the  "hot" 

CO  will  be : 
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where : 


P,  <K) 


Is  Che  return  power  wlchouC  Che  "hoc'* 
CO 


k  • 


Twca 


C*) 


Is  Che  two-way  atmospheric  transmission 
due  to  “hot"  CO 


This  equation  can  be  solved  for  ^4)  for  Che  case  where  /*«  (4)  equals 
Che  system  noise  power: 


T.*,  (.K) 


but : 
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Chen,  Che  maximum  “hot"  CO  number  density  can  be  written: 


if4K 
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where  7^^  is  Che  total  absorption  cross-section  due  to  “hot"  CO. 

In  order  to  calculate  the  maximum  “hot"  CO  density  that  can  be 
detected,  the  total  absorption  cross-section  for  the  “hot"  CO  must  be 
calculated.  Because  of  the  manner  In  which  the  Lldar  equation  was  writ 
ten,  two  separate  cross-sections  contribute  to  this  total;  the  cross- 
section  due  to  Che  R(2)  line  of  CO,  and  Che  cross-section  due  to  all 
other  nearby  CO  lines.  However,  Che  cross -section  due  to  the  R(2)  line 
will  dominate  by  at  least  two  orders  of  magnitude,  and  the  other  term 
can  be  neglected.  Then,  equations  (26),  (46),  and  (47)  can  be  used  to 


calculate  the  "hot"  CO  cross-section,  which  Is  tabulated  as  a  function 


of  altitude  and  temperature  In  Table  VI. 

The  maximum  number  density  that  can  be  detected  can  now  be  calcu¬ 
lated  from  equation  (122).  These  results  are  also  given  In  Table  VI. 

It  should  be  remembered  that  the  calculated  maximum  number  density 
Is  In  addition  to  the  existing  CO  number  density  at  the  altitude  of 
Interest.  For  this  reason.  It  makes  little  sense  to  Investigate  the 
minimum  number  density,  since  the  standard  CO  content  Is  already  well 
above  the  system  threshold. 


Table  VI 

Heterodyne  Detection  of  Hot  CO 


400 

7.949 

X 

10-1« 

3.040  X 

10^^ 

500 

1.544 

X 

10-^7 

1.565  X 

10^2 

600 

2.381 

X 

lo"^^ 

1.015  X 

10^^ 

700 

3.148 

X 

lO"^^ 

7.677  X 

10^^ 

800 

3.882 

X 

10-'^ 

6.226  X 

lO^' 

900 

4.585 

X 

10-'^ 

5.271  X 

io“ 

1000 

5.147 

X 

10-'^ 

4.696  X 

10^1 

400 

8.943 

X 

io-'« 

2.888  X 

10^^ 

500 

1.725 

X 

lo"^^ 

1.497  X 

10^^ 

600 

2.631 

X 

lo”^^ 

9.816  X 

10^^ 

700 

3.578 

X 

10-'^ 

7.218  X 

10^^ 

800 

4.342 

X 

lo'^^ 

5.948  X 

10^^ 

900 

5.038 

X 

10-'^ 

5.126  X 

10^' 

1000 

5.677 

X 

10-^' 

4,549  X 

io“ 

400 

1.001 

X 

lo”^^ 

2.426  X 

10^^ 

500 

1.951 

X 

10“  ^  7 

1.245  X 

10^^ 

600 

2.936 

X 

10-^7 

8.270  X 

10^^ 

700 

3.931 

X 

10-'^ 

6.177  X 

10^^ 

800 

4.915 

X 

10-'^ 

4.940  X 

10^^ 

900 

5.733 

X 

10-'^ 

4.235  X 

10^^ 

1000 

6.317 

X 

10-^' 

3.844  X 

10^' 

VIII.  Conclusions  and  Recommendations 


Conclusions 

On  the  basis  of  the  analysis  presented  In  this  report,  the 
following  conclusions  are  made: 


1.  A  range-resolved  DIAL  Ildar  system  on  board  the  Space 
Shuttle  can  be  used  to  measure  CO  concentrations  in  the 
troposphere,  but  stratospheric  measurements  will  be 
marginal  due  to  large  measurement  errors. 

2.  Heterodyne  detection  of  the  return  signal  is  required. 
Direct  detection  sensitivity  was  found  to  be  unsatisfac¬ 
tory  due  to  thermal  noise  In  the  detection  system. 

3.  Of  the  major  Interfering  species  in  the  atmosphere  at 
4.64  (carbon  dioxide,  nitrous  oxide,  ozone,  and  water- 
vapor),  only  water-vapor  provides  a  significant  contri¬ 
bution  to  the  extinction  coefficient. 

4.  Addition  of  hot  CO  near  the  ground  can  be  detected  up  to  a 
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maximum  concentration  of  about  10  cm  depending  on  the 
altitude  and  temperature. 


RecommendaClons 


Based  on  the  results  of  this  thesis  the  following  recommendations 
are  made : 

1.  A  follow-on  study  be  conducted  to  determine  the  sen¬ 
sitivity  of  the  results  presented  herein  to  the  atmosphere 
model  used.  In  particular,  the  system  sensitivity  to 
changes  In  the  water-vapor  profile  should  be  addressed, 
since  the  system  will  be  required  to  make  measurements 
over  large  ocean  areas. 

2.  A  follow-on  study  concerning  detection  of  hot  gases  near 
the  earth's  surface  should  be  conducted.  This  study  and 
the  one  by  Wank  both  conclude  that  a  DIAL  system  In  orbit 
can  provide  a  detection  capability  for  such  gases.  It 
would  appear  that  such  a  system  could  be  quite  effective 
In  locating  and  characterizing  emissions  from  a  large 
Industrial  complex  or  the  exhaust  plumes  of  gas  turbine 
or  rocket  engines. 
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CO  Line  Positions  and  Line  Strengths 


Rot.  ID 


(ca“^) 


S,„  (cm2  mol~l  cm“l) 


P(24) 


2041.6668 


1.618  X  10-20 


P(23) 


2046.2764 


2.429  X 


P(22) 

P(21) 

P(20) 


2050.8544 


2055.4007 


2059.9150 


3.570  X  10-20 
5.140  X  10-20 
7.248  X  10-19 


P(19) 

P(18) 

P(17) 

P(16) 

P(15) 

P<14) 

P(13) 

P(12) 

P(ll) 

P(10) 


2064.3973 

2068.8473 


2073.2650 


2077.6501 

2082.0026 


2086.3223 


2090.6091 


2094.8627 

2099.0831 


2103.2701 


1.001  X  10-19 

1.352  X  10-19 
1.787  X  10-19 
2.311  X  10-19 
2.922  X  10-19 
3.609  X  10-19 

4.353  X  10-19 
5.123  X  10-19 
5.877  X  10-19 
6.561  X  10-19 


2107.4236 


7.118  X 


2111.5434 


7.486  X 


2115.6294 


7.607  X 


2119.6814 


7.431  X 


2123.6992 


6.927  X 


2127.6828 


6.084  X 


,*•  m*"-*"*  •*^i_^**i**J*^ ,  •  .*•  ,*•  ***  »*■ 


P(3) 

2131.6320 

4.917  X 

10-19 

P(2) 

2135.5466 

3.467  X 

10-19 

P(l) 

2139.4265 

1.800  X 

10-19 

Band  Center 

2143.2715 

- 

R(0) 

2147.0815 

1.834  X 

10-19 

R(l) 

2150.8564 

3.599  X 

10-19 

R(2) 

2154.5960 

5.201  X 

10-19 

R(3) 

2158.3001 

6.557  X 

10-19 

R(4) 

2161.9686 

7.606  X 

10“19 

R(5) 

2165.6014 

8.313  X 

10-19 

R(6) 

2169.1983 

8.670  X 

10-19 

R(7) 

2172.7592 

8.693  X 

10-19 

R(8) 

2176.2839 

8.422  X 

10-19 

R(9) 

2179.7723 

7.909  X 

10-19 

R(10) 

2183.2242 

7.218  X 

10-19 

R(ll) 

2186.6394 

6.411  X 

10-19 

R(12) 

2190.0179 

5.550  X 

10-19 

R(13) 

2193.3595 

4.688  X 

10-19 

R(14) 

2196.6641 

3.867  X 

10-19 

R(15) 

2199.9314 

3.116  X 

10-19 

R(16) 

2203.1614 

2.456  X 

10-19 

R(17) 

2206.3539 

1.892  X 

10-19 

R(18) 

2209.5087 

1.427  X 

10-19 

R(19) 

2212.6257 

1.053  X 

10-19 

R(20) 

2215.7048 

7.610  X 

o 

1 

o 

.  -  .  , 

1  ■ 

101 
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Atmospheric  Parameters 


NTOT(cm“’) 


1.0000 

2.547 

X 

1019 

1.885 

X 

10l2 

0.8870 

2.311 

X 

1019 

1.710 

X 

10l2 

0.7846 

2.093 

X 

10l9 

1.549 

X 

1012 

0.6920 

1.891 

X 

1019 

1.399 

X 

1012 

0.6085 

1.704 

X 

1019 

1.261 

X 

1012 

0.5334 

1.531 

X 

ioi9 

1.133 

X 

1012 

0.4660 

1.373 

X 

1019 

1.016 

X 

1012 

0.4057 

1.227 

X 

1019 

9.078 

X 

lOll 

0.3519 

1.093 

X 

10l9 

8.087 

X 

lOll 

0.3040 

9.711 

X 

10l8 

7.185 

X 

lOll 

0.2615 

8.598 

X 

1018 

6.362 

X 

lOll 

0.2240 

7.585 

X 

10l8 

5.612 

X 

loll 

0.1915 

6.486 

X 

10l8 

4.799 

X 

loll 

0.1636 

5.543 

X 

10l8 

4.101 

X 

loll 

0.1399 

4.738 

X 

1018 

3.506 

X 

lOll 

0.1195 

4.049 

X 

1018 

2.996 

X 

loll 

0.1022 

3.461 

X 

1018 

2.561 

X 

lOll 

0.08734 

2.959 

X 

10l8 

2.189 

X 

loll 

0.07466 

2.529 

X 

10l8 

1.871 

X 

loll 

0.06383 

2. 162 

X 

10l8 

1.600 

X 

lOll 

0.05457 

1.849 

X 

10l8 

1.368 

X 

loll 

0.04667 

1.574 

X 

1018 

1.165 

X 

loll 

0.03995 

1.341 

X 

10l8 

9.922 

X 

lolo 

0.03422 

1.144 

X 

10l8 

8.464 

X 

loio 

0.02933 

9.759 

X 

10l7 

7.221 

X 

lOlO 

0.02516 

8.334 

X 

10l7 

6.166 

X 

lOlO 

0.02160 

7.123 

X 

10l7 

5.270 

X 

lOlO 

0.01855 

6.092 

X 

10l7 

4.507 

X 

lOlO 

0.01595 

5.214 

X 

1017 

3.858 

X 

lOlO 

0.01372 

4. 466 

X 

1017 

3.304 

X 

lOlO 

0.01181 

3.828 

X 

1017 

2.832 

X 

lOlO 
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Aerosol  Extinction  Cefflclents  for  Clear  Atmosphere 

aerosol  ^aerosol 

absorption  scattering 


1.40 

X 

10~2 

1.44 

X 

10-2 

2.84 

X 

10-2 

9.31 

X 

10-3 

9.54 

X 

10-3 

1.89 

X 

10-2 

4.05 

X 

10-3 

4.17 

X 

10-3 

8.22 

X 

10-3 

1.73 

X 

10-3 

1.78 

X 

10-3 

3.51 

X 

10-3 

8.14 

X 

10-^ 

8.38 

X 

10-4 

1.65 

X 

10-3 

5.13 

X 

10-^ 

5.29 

X 

10-4 

1.04 

X 

10-3 

3.74 

X 

10-^ 

3.85 

X 

10-4 

7.59 

X 

10-4 

3.03 

X 

10-^ 

3.12 

X 

10-4 

6.15 

X 

10-4 

2.97 

X 

10-^ 

3.06 

X 

10-4 

6.03 

X 

10-4 

2.95 

X 

10-^ 

3.04 

X 

10-4 

5.99 

X 

10-4 

2.85 

X 

10-^ 

2.93 

X 

10-4 

5.78 

X 

10-4 

2.72 

X 

10-^ 

2.81 

X 

10-4 

5.53 

X 

10-4 

2.70 

X 

10-^ 

2.78 

X 

10-4 

5.48 

X 

10-4 

2.66 

X 

10-^ 

2.74 

X 

10-4 

5.40 

X 

10-4 

2.54 

X 

10-4 

2.61 

X 

10-4 

5.15 

X 

10-4 

2.43 

X 

10-4 

2.50 

X 

10-4 

4.93 

X 

10-4 

2.30 

X 

10-4 

2.36 

X 

10-4 

4.66 

X 

10-4 

2.23 

X 

10-4 

2.29 

X 

10-4 

4.52 

X 

10-4 

2.17 

X 

10-4 

2.24 

X 

10-4 

4.41 

X 

10-4 

1.97 

X 

10-4 

2.02 

X 

10-4 

3.99 

X 

10-4 

1.54 

X 

10-4 

1.59 

X 

10-4 

3.13 

X 

10-4 

1.13 

X 

10"^ 

1.16 

X 

10-^ 

2.29 

X 

10-^ 

8.32 

X 

10-5 

8.57 

X 

10-5 

1.69 

X 

10-^ 

6.31 

X 

10-5 

6.50 

X 

10-5 

1.28 

X 

10-^ 

4.92 

X 

10-5 

5.07 

X 

10-5 

9.99 

X 

10-5 

4.02 

X 

10-5 

4.14 

X 

10-5 

8.16 

X 

10-5 

2.02 

X 

10-5 

2.08 

X 

10-5 

4.10 

X 

10-5 

2.02 

X 

10-5 

2.08 

X 

10-5 

4.10 

X 

10-5 

2.02 

X 

10-5 

2.08 

X 

10-5 

4.10 

X 

10-5 

2.02 

X 

10-5 

2.08 

X 

10-5 

4.10 

X 

10-5 

2.02 

X 

10-5 

2.08 

X 

10-5 

4.10 

X 

10-5 
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Rayleigh  Atmospheric  Extinction  Coefficient 


(km“^) 

h  (km) 

(km-1) 

2.290 

X 

10-6 

21 

1.415  X 

10-7 

2.078 

X 

10-6 

22 

1.206  X 

10-7 

1.882 

X 

10-6 

23 

1.029  X 

10-7 

1.700 

X 

10-6 

24 

8.776  X 

10-8 

1.532 

X 

10-6 

25 

7.494  X 

10-8 

1.377 

X 

10-6 

26 

6.405  X 

10-8 

1.235 

X 

10-6 

27 

5.478  X 

10-8 

1.103 

X 

10-6 

28 

4.689  X 

10-8 

9.829 

X 

10-7 

29 

4.016  X 

10-8 

8.733 

X 

10-7 

30 

3.422  X 

10-8 

7.732 

X 

10-7 

6.821 

X 

10-7 

5.833 

X 

10-7 

4.985 

X 

10“7 

4.261 

X 

10-7 

3.641 

X 

10-7 

3.112 

X 

10-7 

2.661 

X 

10-7 

2.274 

X 

10-7 

1.944 

X 

10-7 

1.663 

X 

10-7 

h  (km) 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
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1.689 

X 

10-8 

9.23 

X 

10-6 

9.247 

X 

10-6 

1.400 

X 

10-8 

6.82 

X 

10-6 

6.834 

X 

10-6 

1.228 

X 

10-8 

5.17 

X 

10-6 

5.182 

X 

10-8 

1.048 

X 

10-8 

4.03 

X 

10-6 

4.040 

X 

10-6 

8.945 

X 

10-9 

3.29 

X 

10-6 

3.299 

X 

10-8 

7.645 

X 

10-9 

1.66 

X 

10-6 

1.668 

X 

10-8 

6.539 

X 

10-9 

1.66 

X 

10-6 

1.667 

X 

10-8 

5.597 

X 

10-9 

1.66 

X 

10-6 

1.666 

X 

10'8 

4.794 

X 

10-9 

1.66 

X 

10-6 

1.665 

X 

10-8 

4.109 

X 

10-9 

1.66 

X 

10-6 

1.664 

X 

10-8 

vv- 

%*  V 


Appendix  G 

Absorption  Cross-Section  and  Atmospheric  Extinction  Coefficient  of  CO 


(k.m~ 


2.595  X 

10~21 

4.892 

X 

10-4 

1.888  X 

10-21 

3.228 

X 

10-4 

1.513  X 

10-21 

2.344 

X 

10-4 

1.203  X 

10-21 

1.683 

X 

10-4 

9.478  X 

10-22 

1.195 

X 

10-4 

7.395  X 

10-22 

8.379 

X 

10-5 

5.710  X 

10-22 

5.801 

X 

10-5 

4.971  X 

10-22 

4.513 

X 

10-5 

3.289  X 

10-22 

2.660 

X 

10-5 

2.449  X 

10-22 

1.760 

X 

10-5 

0 

1.803  X 

10-22 

1.147 

X 

10-5 

1 

1.304  X 

10-22 

7.318 

X 

10-6 

2 

1.112  X 

10-22 

5.336 

X 

10-6 

3 

9.497  X 

10-23 

3.895 

X 

10-6 

4 

8.121  X 

10-23 

2.847 

X 

10-6 

5 

6.937  X 

10-23 

2.078 

X 

10-6 

6 

5.933  X 

10-23 

1.519 

X 

10-6 

7 

5.070  X 

10-23 

1.110 

X 

10-6 

8 

4.334  X 

10-23 

8.109 

X 

10-7 

9 

3.705  X 

10-23 

5.928 

X 

10-7 

;o 

3.168  X 

10-23 

4.334 

X 

10-7 

2.775  X 

10-23 

3.233  X 

10-7 

2.440  X 

10-23 

2.421  X 

10-7 

2.145  X 

10-23 

1.816  X 

10-7 

1.887  X 

10-23 

1.363  X 

10-7 

1.657  X 

10-23 

1.022  X 

10-7 

1.459  X 

10-23 

7.689  X 

10-8 

1.285  X 

10-23 

5.791  X 

10-8 

1.133  X 

10-23 

4.371  X 

10-8 

9.990  X 

10-24 

3.301  X 

10-8 

8.812  X 

10-24 

2.496  X 

10-8 
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Absorption  Cross -Section  of  Water-Vapor  at  Sea-Level 


V,  (cm~l) 


(cm 


~lmol“lc 


2141.508 

0.235  X 

o 

N> 

1.395 

X 

10-30 

2141.519 

0. 145  X 

10“23 

8.625 

X 

10-29 

2142.246 

0.247  X 

10-23 

1.647 

X 

10-28 

2142.301 

0. 336  X 

10-24 

2.261 

X 

10-29 

2142.495 

0.228  X 

10-25 

1.584 

X 

10-30 

2143.369 

0.129  X 

10-25 

1.041 

X 

10-30 

2143.982 

0.162  X 

10-25 

1.462 

X 

10-30 

2144.708 

0.101  X 

10-23 

1.050 

X 

10-28 

2144.808 

0.742  X 

10-23 

7.874 

X 

10-28 

2145.468 

0.579  X 

10-22 

7.064 

X 

10-27 

2145.639 

0.366  X 

10-25 

4.637 

X 

10-30 

2145.678 

0.534  X 

10-23 

6.825 

X 

10-28 

2146.727 

0.157  X 

10-25 

2.577 

X 

10-30 

2147.337 

0.215  X 

10-25 

4.146 

X 

10-30 

2147.406 

0.304  X 

10-22 

5.975 

X 

10-27 

2148.187 

0.368  X 

10-23 

9.100 

X 

10-28 

2148.342 

0.448  X 

10-23 

1.163 

X 

10-27 

2148.479 

0.208  X 

10-24 

5.645 

X 

10-29 

2148.831 

0.130  X 

10-24 

3.972 

X 

10-29 

2148.962 

0.726  X 

10-25 

2.322 

X 

10-29 

2149.130 

0.103  X 

10-23 

3.500 

X 

10-28 

2149.620 

0.561  X 

10-25 

2.299 

X 

10-29 

2151.194 


2151.423 


2151.528 


2151.529 


2152.559 


2152.773 


0.563  X  10-23 
0.438  X  10-25 
0.219  X  10-25 
0.656  X  10-25 
0.451  X  10-22 
0.503  X  10-25 


4.928  X  10-27 
4.406  X  10-29 
2.356  X  10-29 
7.061  X  10-29 
1.097  X  10-25 
1.526  X  10-28 


*Cs^*w*’** 


'K'fy. 


2152.861 


0.168  X  10-25 


5.624  X  10-29 


2152.922 


2153.241 


0.113  X  10-25 
0.189  X  10-25 


4.062  X  10-29 
1.034  X  10-28 


2153.537 


0.366  X  10-25 


3.265  X  10-28 


2154.056 


2154.355 


2154.710 


2155.000 


0.981  X  10-25 
0.299  X  10-24 
0.488  X  10-23 
0.101  X  10-23 


3.304  X  10-27 
4.794  X  10-26 
4.125  X  10-24 
6.551  X  10-26 


•A  . 
•  .V 

«..  fc." 


2155.217 


2155.295 


2156.565 


2156.636 


0.152  X  10-25 
0.908  X  10-25 
0.107  X  10-22 
0.582  X  10-2^ 


4.122  X  10-28 
1.938  X  10-27 
2.836  X  10-26 
1.437  X  10-27 


»  •  - 
•  I 

••  -V-'. 


2157.454 


0.108  X  10-24 


1.355  X  10-28 


2158.105 


0.127  X  10-23 


1.056  X  10-27 


2159.842 


0.170  X  10-25 


6.313  X  10-30 


2160.006 


0.165  X 


5.761  X  10-28 


2160.263 


0.557  X  10-25 


1.772  X  10-29 


2160.697 


0.429  X 


1. 177  X 


2161.427 


2161.726 


0.140  X 
0.253  X 


3.064  X  10"30 
5.082  X  10-27 


t: 


2161.727 

0.759  X 

10-22 

1.524 

X 

10-26 

2162.012 

0.195  X 

10"23 

3.620 

X 

10-28 

2162.147 

0.205  X 

10-25 

3.671 

X 

10-30 

2162.891 

0.810  X 

10-23 

1.202 

X 

10-27 

2163.430 

0.242  X 

10-22 

3.165 

X 

10-27 

2165.223 

0.110  X 

10-23 

9.938 

X 

10-29 

2166.112 

0.472  X 

10-24 

3.631 

X 

10-29 

2166.380 

0.117  X 

10-25 

8.595 

X 

10-31 

2167.346 

0.378  X 

10-23 

2.372 

X 

10-28 

2167.771 

0.310  X 

10-24 

1.822 

X 

10-29 

TOTAL  -  4.435  x  10“24 


[•jF 


Atmospheric  Extinction  Coefficient  of  Water-Vapor 


h  (km) 

%20  (cm“3) 

(cm2) 

,  (km-1) 

0 

1.919  X 

10l7 

4.435 

X 

10-24 

8.511 

X 

10-2 

1 

1.375  X 

10l7 

2.703 

X 

10-24 

3.717 

X 

10-2 

2 

9.567  X 

10l6 

1.988 

X 

10-24 

1.902 

X 

10-2 

3 

* 

6.521  X 

ioi6 

1.443 

X 

10-24 

9.410 

X 

10-3 

A 

3.474  X 

10l6 

1.034 

X 

10-24 

3.592 

X 

10-3 

5 

* 

2.349  X 

10l6 

7.297 

X 

10-25 

1.714 

X 

10-3 

6 

1.223  X 

10l6 

5.071 

X 

10-25 

6.202 

X 

10-^ 

7 

* 

8.013  X 

1015 

3.466 

X 

10-25 

2.777 

X 

10-4 

8 

3.796  X 

10l5 

2,326 

X 

10-25 

8.829 

X 

10-5 

9 

* 

2.197  X 

10l5 

1.531 

X 

10-25 

3.364 

X 

10-5 

10 

5.972  X 

lOl'^ 

9.905 

X 

10-26 

5.915 

X 

10-6 

11 

* 

3.575  X 

lOl^ 

6.237 

X 

10-26 

2.230 

X 

10-6 

12 

1.178  X 

loi^ 

5.306 

X 

10-26 

6.250 

X 

10-7 

13 

* 

7.147  X 

10l3 

4.533 

X 

10-26 

3.240 

X 

10-7 

14 

2.514  X 

10l3 

3.877 

X 

10-26 

9.747 

X 

10-8 

15 

* 

2.175  X 

10l3 

3.311 

X 

10-26 

7.201 

X 

10-8 

16 

1.836  X 

10l3 

2.832 

X 

10-26 

5.200 

X 

10-8 

17 

* 

1.589  X 

10l3 

2.420 

X 

10-26 

3.845 

X 

10-8 

18 

1.342  X 

10l3 

2.069 

X 

10-26 

2.777 

X 

10-8 

19 

* 

1.340  X 

10l5 

1.769 

X 

10-26 

2.370 

X 

10-8 

20 

1.338  X 

10l3 

1.512 

X 

10-26 

2.023 

X 

10-8 

*  Indicates 

a 

linearly  Interpolated 

value  of 

number  density 

V  r. ■ 


V 


. 
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21 


22 

23 

24 

25 

26 

27 

28 

29 

30 


*  1.445  X  10l3 

1.552  X  10l3 

*  1.686  X  10l3 

1.820  X  10l3 

*  1.975  X  10l3 

2.130  X  10^3 

*  1.828  X  10l3 

1.526  X  10l3 

*  1.305  X  10l3 

1.083  X  10l3 


1.351 

X 

10-26 

1.952 

X 

10-8 

1.214 

X 

10-26 

1.884 

X 

10-8 

1.091 

X 

10-26 

1.839 

X 

10-8 

9.799 

X 

10-27 

1.783 

X 

10-8 

8.768 

X 

10-27 

1.732 

X 

10-8 

7.885 

X 

10-27 

1.680 

X 

10-8 

7.091 

X 

10-27 

1.296 

X 

10-8 

6.38i 

X 

10-27 

9.737 

X 

10-9 

5.742 

X 

10-27 

7.493 

X 

10-9 

5.169 

X 

10-27 

5.598 

X 

10-9 

*  Indicates  a  linearly  Interpolated  value  of  number  density 


Appendix  K 

Absorption  Cross-Section  of  CO  R(2)  Line  at  2154.6050  cm  ^ 


(km) 

Sm  (cm“^mol“^cm2) 

o<  (cffl-l) 

Vco  (cm2) 

0 

5,201  X  10-19 

7.0  X  10-2 

2.327  X  10-18 

1 

4.173  X  10-19 

6.4  X  10-2 

2.035  X  10-18 

2 

3.744  X  10-19 

5.7  X  10-2 

2.040  X  10-18 

3 

3.339  X  10-19 

5.1  X  10-2 

2.021  X  10-18 

4 

2.960  X  10-19 

4.5  X  10-2 

2.013  X  10-18 

5 

2.606  X  10-19 

4.0  X  10-2 

1.974  X  10-18 

6 

2.277  X  10-19 

3.6  X  10-2 

1.895  X  10-18 

7 

1.975  X  10-19 

3.1  X  10-2 

1.870  X  10-18 

8 

1.697  X  10-19 

2.8  X  10-2 

1.749  X  10-18 

9 

1.446  X  10-19 

2.4  X  10-2 

1.681  X  10-18 

10 

1.222  X  10-19 

2.1  X  10-2 

1.565  X  10-18 

11 

1.019  X  10-19 

1.8  X  10-2 

1.442  X  10-18 

12 

1.016  X  10-19 

1.6  X  10-2 

1.535  X  10-18 

13 

1.016  X  10-19 

1.3  X  10-2 

1.682  X  10-18 

14 

1.016  X  10-19 

1.1  X  10-2 

1.761  X  10-18 

15 

1.016  X  10-19 

9.8  X  10-3 

1.790  X  10-18 

16 

1.016  X  10-19 

8.4  X  10-3 

1.792  X  10-18 

17 

1.016  X  10-19 

7.1  X  10-3 

1.747  X  10-18 

18 

1.016  X  10-19 

6.1  X  10-3 

1.669  X  10-18 

19 

1.016  X  10-19 

5.2  X  10-3 

1.557  X  10-18 

20 

1.016  X  10-19 

4.5  X  10-3 

1.437  X  10-18 

1.043 

X 

10-19 

3.8 

X 

10-3 

1.322 

X 

10-18 

1.073 

X 

10-19 

3.3 

X 

10-3 

1.227 

X 

10-18 

1.104 

X 

10-19 

2.8 

X 

10-3 

1.108 

X 

10-18 

1.135 

X 

10-19 

2.4 

X 

10-3 

9.994 

X 

10-19 

1.163 

X 

10-19 

2.0 

X 

10-3 

8.710 

X 

10-19 

1.196 

X 

10-19 

1.7 

X 

10-3 

7.715 

X 

10-19 

1.229 

X 

10-19 

1.5 

X 

10-3 

7.049 

X 

10-19 

1.262 

X 

10-19 

1.3 

X 

10-3 

6.315 

X 

10-19 

1.296 

X 

10-19 

1.1 

X 

10-3 

5.520 

X 

10-19 

1.331 

X 

10-19 

9.5 

X 

10-^ 

4.914 

X 

10-19 

Appendix  L 


Atmospheric  Transmission  at  2154.6050  cm~^  IXie  to  R(2)  Line  of  CO 


(km) 

t/.  (R) 

1, 

0 

0.4159 

0.001847 

1 

0.4986 

0.004441 

2 

0.5315 

0.008907 

3 

0.5681 

0.01676 

4 

0.6019 

0.02950 

5 

0.6393 

0.04901 

6 

0.6804 

0.07666 

7 

0.7121 

0.1127 

8 

0.7536 

0.1582 

9 

0.7854 

0.2099 

10 

0.8194 

0.2673 

11 

0.8527 

0.3262 

12 

0.8630 

0.3826 

13 

0.8711 

0.4433 

14 

0.8838 

0. 5089 

15 

0.8983 

0.5758 

16 

0.9123 

0.6410 

17 

0.9264 

0. 7026 

18 

0.9395 

0. 7585 

19 

0.9514 

0.8073 

20 

0.9614 

0.8485 

i 

:.‘>N 


V  V 


I'.'" 


6. 138  X 

10-19 

0.9858 

0.9169 

6.390  X 

10-19 

0.9874 

0.9301 

6.594  X 

10-19 

0.9889 

0.9420 

6.819  X 

10-19 

0.9902 

0.9526 

7.004  X 

10-19 

0.9914 

0.9620 

7.238  X 

10-19 

0.9924 

0.9704 

7.458  X 

10-19 

0.9933 

0.9778 

7.669  X 

10-19 

0.9941 

0.9844 

7.890  X 

10-19 

0.9948 

0.9902 

8.140  X 

10-19 

0.9954 

0.9954 
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